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Abstract 


Studies of the spin, parity and tensor couplings of the Higgs boson in the H ZZ* Ai, 
H WW* evfiv and H ^ yy decay processes at the LHC are presented. The invest¬ 
igations are based on 25 fb“^ of pp collision data collected by the ATLAS experiment at 
= 7 TeV and ^fs = % TeV. The Standard Model (SM) Higgs boson hypothesis, corres¬ 
ponding to the quantum numbers = 0"^, is tested against several alternative spin scenarios, 
including non-SM spin-0 and spin-2 models with universal and non-universal couplings to 
fermions and vector bosons. All tested alternative models are excluded in favour of the SM 
Higgs boson hypothesis at more than 99.9% confidence level. Using the H ZZ* —> At and 
H WW* ^ evpv decays, the tensor structure of the interaction between the spin-0 boson 
and the SM vector bosons is also investigated. The observed distributions of variables sens¬ 
itive to the non-SM tensor couplings are compatible with the SM predictions and constraints 
on the non-SM couplings are derived. 
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1 Introduction 


The discovery of a Higgs boson by the ATLAS [1] and CMS [2] experiments at the Large Hadron Collider 
(LHC) at CERN marked the beginning of a new era of experimental studies of the properties of this 
new particle. In the Standard Model (SM), the Higgs boson is a CP-even scalar particle, = O'*"''. ^ 
Theories of physics beyond the SM (BSM) often require an extended Higgs sector featuring several 
neutral Higgs bosons. Such cases may include CP-mixing in the Higgs boson interactions, which could 
result in observable differences in the kinematics of final-state particles produced in their decays. A 
review of the phenomenology in the determination of Higgs boson spin and CP properties can be found 
in Ref. [3] and references therein. 

Previous determinations of the Higgs boson spin and CP quantum numbers by the ATLAS and CMS 
Collaborations are reported in Refs. [4] and [5]. Results on the same subject have also been published by 
the DO and CDF Collaborations in Ref. [6]. All these studies indicate the compatibility of the spin and 
CP properties of the observed Higgs boson with the SM predictions. The ATLAS measurement excluded 
several alternative spin and parity hypotheses in favour of the quantum numbers predicted by the SM. In 
addition to the exclusion of several non-SM spin hypotheses, the CMS measurement probed the tensor 
structure of the Higgs boson decay to SM vector bosons in the spin-0 scenario. This paper complements 
the previous ATLAS study of the Higgs boson spin and parity. The new study takes advantage of im¬ 
provements to the analysis strategy and to the modelling used to describe alternative spin hypotheses, and 
includes studies on CP-mixing for the spin-0 scenario. The improved theoretical framework is based on 
the Higgs boson characterisation model described in Refs. [3, 7]. 

The study of the spin and parity properties of the Higgs boson presented in this paper is based on the 
H yy, H ZZ* ^ 4i and H —> WW* evfiv decay channels and their combination. The 
H —> WW* evjjv analysis is described in detail in a separate publication [8]. These analyses are based 
on 4.5 fb“' and 20.3 fb“^ of pp collision data collected by the ATLAS experiment at centre-of-mass 
energies of 7 TeV and 8 TeV, respectively. For the H WW* —> evpv studies only the data collected 
at a centre-of-mass energy of 8 TeV are used. The SM hypothesis - O'*” is compared to alternative 
spin-0 models: a pseudoscalar boson = 0~ and a BSM scalar boson - Oj^ [9, 10], which describes 
the interaction of the Higgs boson with the SM vector bosons with higher-dimension operators discussed 
in Section 3.1. Graviton-like tensor models with - 2* with universal and non-universal couplings [3, 
7] are also considered. In these tests of fixed spin and parity hypotheses it is assumed that the resonance 
decay involves only one CP eigenstate. 

In addition to the fixed spin and parify hypofhesis fesfs, fhe possible presence of BSM terms in fhe Fag- 
rangian describing fhe HVV vertex^ of fhe spin-0 resonance is also invesfigafed. The HVV inferacfion 
is described in terms of an effecfive Fagrangian fhaf confains fhe SM inferacfion and BSM CP-odd and 
CP-even terms [3, 7]. The relative fractions of the CP-odd and CP-even BSM contributions to the ob¬ 
served Higgs boson decays are constrained, and limits on the corresponding BSM tensor couplings are 
derived. 

This paper is organised as follows. In Section 2 the ATFAS detector is described. In Section 3 the 
theoretical framework used to derive the spin and parity models, as well as the parameterisation used 
to describe the HVV coupling tensor structure, are discussed. In Section 4, the choice of Monte Carlo 
generators for the simulation of signal and backgrounds is described. The analyses of fixed spin and parify 

* In the following, for brevity, only the J’’ label is used to indicate the spin and CP quantum numbers. 

^ In this paper the symbol V is used to describe a massive SM vector boson, namely either a IV or a Z boson. 
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hypotheses for the three decay channels and their combination are presented in Section 5. Individual and 
combined studies of the tensor structure of the HVV interaction are presented in Section 6 . Concluding 
remarks are given in Section 7. 


2 The ATLAS detector 

The ATLAS detector is described in detail in Ref. [11]. ATLAS is a multi-purpose detector with a forward- 
backward symmetric cylindrical geometry. It uses a right-handed coordinate system with its origin at the 
nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The 
x-axis points from the IP to the centre of the LHC ring, and the i/-axis points upward. Cylindrical co¬ 
ordinates (r, 4>) are used in the transverse plane, ^ being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined as 77 = - lntan(0/2), where G is the polar angle. 

At small radii from the beamline, the inner detector (ID), immersed in a 2 T magnetic field produced by 
a thin superconducting solenoid located in front of the calorimeter, is made up of fine-granularity pixel 
and microstrip detectors. These silicon-based detectors cover the range \t]\ <2.5. A gas-filled straw-tube 
transition-radiation tracker (TRT) complements the silicon tracker at larger radii and also provides elec¬ 
tron identification based on transition radiation. The electromagnetic (EM) calorimeter is a lead/liquid- 
argon sampling calorimeter with an accordion geometry. The EM calorimeter is divided into a barrel 
section covering |? 7 | < 1.475 and two end-cap sections covering 1.375 < |? 7 | < 3.2. Eor |? 7 | < 2.5 it is 
divided into three layers in depth, which are finely segmented in 77 and (j). An additional thin presampler 
layer, covering I 77 I < 1 . 8 , is used to correct for fluctuations in energy losses of particles before they reach 
the calorimeter. Hadronic calorimetry in the region I 77 I <1.7 uses steel absorbers and scintillator tiles 
as the active medium. Eiquid argon with copper absorbers is used in the hadronic end-cap calorimeters, 
which cover the region 1.5 < I 77 I < 3.2. A forward calorimeter using copper or tungsten absorbers with 
liquid argon completes the calorimeter coverage up to I 77 I = 4.9. The muon spectrometer (MS) measures 
the deflection of muon trajectories with I 77 I < 2.7, using three stations of precision drift tubes, with cath¬ 
ode strip chambers in the innermost layer for |? 7 | > 2.0. The deflection is provided by a toroidal magnetic 
field with an integral of approximately 3 Tm and 6 Tm in the central and end-cap regions of the ATEAS 
detector, respectively. The muon spectrometer is also instrumented with dedicated trigger chambers, the 
resistive-plate chambers in the barrel and thin-gap chambers in the end-cap, covering I 77 I < 2.4. 


3 Theoretical models 

In this section, the theoretical framework for the measurements of the spin and parity of the resonance 
is discussed. An effective field theory (EET) approach is adopted to describe the interaction between 
the resonance and the SM vector bosons, following the Higgs boson characterisation model described in 
Refs. [3, 7]. Three possible BSM scenarios for the spin and parity of the boson are considered: 

• the observed resonance is a spin -2 particle, 

• the observed resonance is a pure BSM spin-0 CP-even or CP-odd Higgs boson, 

• the observed resonance is a mixture of the SM spin-0 state and a BSM spin-0 CP-even or CP-odd 
state. 
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The third case would imply CP-violation in the Higgs sector. In the case of CP mixing, the Higgs boson 
would be a mass eigenstate, but not a CP eigenstate. In all cases, only one resonance with a mass of 
about 125 GeV is considered. It is also assumed that the total width of the resonance is small compared 
to the typical experimental resolution of the ATLAS detector (of the order of 1-2 GeV in the four-lepton 
and yy final states, as documented in Ref. [12]). Interference effects between the BSM signals and SM 
backgrounds are neglected. 

The EFT approach, used by the Higgs boson characterisation model, is only valid up to a certain energy 
scale, A. The models described in Ref. [7] assume that the resonance structure corresponds to one new 
boson (X(7^) with = 0* or 2'^), assuming that any other BSM particle only exists at an energy scale 
larger than A. The A scale is set to 1 TeV to account for the experimental results obtained at the LHC and 
previous collider experiments, which do not show any evidence of new physics at lower energy scales. 

The case where the observed resonance has = 1* is not studied in this paper. The H ^ yy decay is 
forbidden by the Landau-Yang theorem [13, 14] for a spin-1 particle. Moreover, the spin-1 hypothesis 
was already studied in the previous ATLAS publication [4] in the H ZZ* —> Ai and H —> WW* —> evjiv 
decays and excluded at a more than 99% confidence level. 


3.1 The spin-0 hypothesis 

In the spin-0 hypothesis, models with fixed spin and parity, and models with mixed SM spin-0 and BSM 
spin-0 CP-even and CP-odd contributions are considered. In Ref. [7], the spin-0 particle interaction with 
pairs of IT or Z bosons is given through the following interaction Lagrangian: 

£q = {cos(a)A-sM \\gHZzZf^Z^^ + 

- — [cos(a)/c-//zzZ^yZ^’' + sin(Qr)/fAzz2^vZ^''] (1) 

[cos(a)^™W;,lT-^" + sm{a)KAwwW;,W-^''\^XQ. 

Here represents the vector-boson field (V = Z, W-), the V^'' are the reduced field tensors and the 
dual tensor is defined as = ^e^^P^'Vpo-. The symbol A denotes the EFT energy scale. The symbols 
khvv and kavv denote the coupling constants corresponding to the interaction of the SM, BSM 
CP-even or BSM CP-odd spin-0 particle, represented by the Xq field, with ZZ or WW pairs. To ensure 
that the Lagrangian terms are Hermitian, these couplings are assumed to be real. The mixing angle a 
allows for production of CP-mixed states and implies CP-violation for a 7 !^ 0 and a + n, provided the 
corresponding coupling constants are non-vanishing. The SM couplings, are proportional to the 

square of the vector boson masses: gnvv ^y. Other higher-order operators described in Ref. [7], 
namely the derivative operators, are not included in Eq. (1) and have been neglected in this analysis since 
they induce modifications of the discriminant variables well below the sensitivity achievable with the 
available data sample. 

As already mentioned, for the spin-0 studies the SM Higgs boson hypothesis is compared to two altern¬ 
atives: the CP-odd - 0~ and the BSM CP-even hypotheses. All three models are obtained 

by selecting the corresponding parts of the Lagrangian described in Eq. (1) while setting all other con¬ 
tributions to zero. The values of the couplings corresponding to the different spin-0 models are listed in 
Table 1. 

The investigation of the tensor structure of the HVV interaction is based on the assumption that the ob¬ 
served particle has spin zero. Following the parameterisation defined in Eq. (1), scenarios are considered 
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7^ 

Model 

Values of tensor couplings 

/rsM Khvv kavv « 

0+ 

SM Higgs boson 

1 

0 

0 

0 

% 

BSM spin-0 CP-even 

0 

1 

0 

0 

0~ 

BSM spin-0 CP-odd 

0 

0 

1 

njl 


Table 1; Parameters of the benchmark scenarios for spin-0 boson tensor couplings used in tests (see Eq. (1)) of the 
fixed spin and parity models. 


where only one CP-odd or one CP-even BSM contribution at a time is present in addition to the SM 
contribution. To quantify the presence of BSM contributions m H ^ ZZ* and H —> WW* decays, the 
ratios of couplings • tana and K-hvvIi^sm are measured. Here kavv and khvv are defined as 

follows: 

^ 1 V 1 V 

T^avv - and khvv - (2) 

4 A 4 A 

where v is the vacuum expectation value [15] of the SM Higgs field. 

The mixing parameters {kavvIksm) • tana and khvvIksm correspond to the ratios of tensor couplings 
g^lgi and g 2 lgi proposed in the anomalous coupling approach described in Refs. [9, 10]. To compare 
the results obtained in this analysis to other existing studies, the final resulfs are also expressed in terms 
of the effective cross-section fractions (,fg 2 ,<Pg 2 ) and {fgA,(t>gA) proposed in Refs. [3] and [9, 10]. Further 
details of these conversions are given in Appendix A. 

The BSM terms described in Eq. (1) are also expected to change the relative contributions of the vector- 
boson fusion (VBF) and vector-boson associated production {VH) processes with respect to the gluon- 
fusion (ggF) production process, which is predicted to be the main production mode for the SM Higgs 
boson at the FHC. For large values of the BSM couplings, at the FHC energies, the VBF production 
mode can have a cross section that is comparable to the ggF process [16]. This study uses only kinematic 
properties of particles from H VV* decays to derive information on the CP nature of the Higgs boson. 
The use of the signal rate information for different production modes, in the context of the EFT analysis, 
may increase the sensitivity to the BSM couplings at the cost of a loss in generality. For example the ratio 
of the VBF and VH production modes with respect to the ggF one can be changed by a large amount for 
non-vanishing values of the BSM couplings. In the studies presented in this paper the predictions of the 
signal rates are not used to constrain the BSM couplings. 

As described in Section 6.2, only events with no reconstructed jets (the 0-jet category) are used in the H 
WW* evjJLV analysis for the studies of the tensor structure; hence this analysis has little sensitivity to 
the VBF production mode. The H ZZ* —> analysis also has little sensitivity to this production mode 

since it is mainly based on variables related to the four-lepton kinematics. The Boosted Decision Tree 
(BDT) algorithm [17] used to discriminate signals from the ZZ* background, described in Sections 5.4 
and 6.3, includes the transverse momentum of the four-lepton system and is trained on simulated samples 
of ggF-produced signals. An enhancement of the VBF production mode would improve the separation 
between background and signal since it predicts larger values of the transverse momentum spectrum for 
events produced via VBF than via ggF [3]. 
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3.2 The spin-2 hypothesis 


In the Higgs boson characterisation model [7], the description of the interaction of a spin-2 particle with 
fermions and vector bosons is described by the following Lagrangian: 



(3) 


The spin-2 tensor field is chosen to interact with the energy-momentum tensors, and T^y, of 
any vector boson V and fermion /, as inspired by gravitation theories. The strength of each interaction 


is determined by the couplings /cy and Kf. In the simplest formulation, all couplings are equal. This 
scenario is referred to as universal couplings (UC), while scenarios with different values of the couplings 
are referred to as non-universal couplings (non-UC). In the UC scenario, the production of a spin-2 
particle in pp collisions is expected to be dominated by QCD processes, with negligible contributions 
from electroweak (EW) processes (i.e. from processes involving EW boson propagators). Simulation 
studies based on MADGRAPH5_aMC@NEO [16] , which implements the Eagrangian described in Eq. (3), 
predict for the production cross section in the UC scenario ctew/ctqcd - 3 x 10~^. These studies also 
show that EW production of the spin-2 resonance would occur mainly in association with a massive 
EW boson {WX, ZX). Present observations do not show a dominant VH production mechanism, hence 
suggesting that ctew is significantly smaller than ctqcd- This paper considers only QCD production for 
all the spin-2 benchmark scenarios. 

The UC models predict a branching ratio of about 5% to photon pairs and negligible branching ratios 
to massive EW gauge boson pairs, WW* and ZZ*. This prediction is disfavoured by the experimental 
measurements [ 1 8-20] and therefore the equality between all couplings k cannot hold. In the benchmark 
scenarios studied in this paper, each of the couplings Ky^/, kz, and Ky is assumed to be independent of 
all the other couplings. In the following, the UC scenario only refers to Kq - Kg, without implying the 
equality for the other k values. 

The simplest QCD production processes, gg ^ X and qq ^ X (where q refers to light quarks), yield 
different polarisations for the spin-2 particle X, and hence different angular distributions of its decay 
products. These mechanisms are considered in the model of a graviton-like tensor with minimal coup¬ 
lings proposed in Refs. [9, 10], which has been studied experimentally in Ref. [4]. The EET Eagrangian, 
however, also allows for more complex processes with emission of one or more additional partons. Eor 
instance, processes with one-parton emission, like qg —> qX and qg —> qX, can produce a spin-2 state 
through either a qqX or a ggX vertex. When two partons are emitted, as in gg qqX or qq qqX, 
the spin-2 production may occur through qqX or ggX vertices, respectively, such that the polarisation of 
X is not uniquely determined by the initial state. Moreover, the EET also allows for four-leg vertices 
like qqgX. These additional diagrams effectively change the polarisation of the particle X, compared to 
what is assumed by the model in Refs. [9, 10]. As a consequence, the angular distributions of the decay 
products become harder to separate from those expected for a scalar resonance. 

The QCD production of a spin-2 particle is driven by the values of the couplings Kg, Kq. Presently, there 
are no experimental constraints on the ratio Kg! Kg from observed decay modes, since the separation of jets 
initiated by gluons or by light quarks is experimentally difficult and has not yet been attempted in Higgs 
boson studies. The ratio KqlKg can thus be regarded as a free parameter. When Kg t Kg, the spin-2 model 
predicts an enhancement of the tail of the distribution of the transverse momentum, p^, of the spin-2 
particle. Such a high-p^ tail is not present for the Kg = Kg (UC) case. As stated before, however, the EETs 


6 




are valid only up to some energy scale, A. At higher energies, new physics phenomena are expected to 
enter to regularise the anomalous ultra-violet behaviour. 

In the present analysis, a selection pj < 300 GeV is applied when investigating non-UC scenarios, 
Kq t K-q. In addition, for the non-UC scenarios, analyses using a tighter selection < 125 GeV are also 
performed. This is a conservative choice for the selection, as the EFT must describe the physics at 
least up to the mass of the observed resonance. It has been verified that the choice of the Pj selection does 
not affect the results for the UC scenario. Even assuming the Pj < 300 GeV selection, some choices of 
Kq!Kg produce high-p^ tails incompatible with the observed differential distribution reported in Refs. [21, 
22]. For this reason the investigated range of the KqjKg ratio is limited to between zero and two. The spin-2 
scenarios considered in this study are presented in Table 2. The Kq - Kg model is referred to hereafter as 
the UC scenario. The Kq = 0 case implies a negligible coupling to light quarks, whereas the Kq = iKg case 
is an alternative scenario with an enhanced coupling to quarks. 


Values of spin-2 quark and gluon couplings 

p^ selections (GeV) 

II 

Universal couplings 

- 

- 

Kq = ^ 

Eow light-quark fraction 

< 300 

< 125 

Kq — '^Kg 

Eow gluon fraction 

< 300 

< 125 


Table 2: Choices of the couplings to quarks Kq and to gluons Kg studied for the spin-2 benchmark scenarios. The 
values of the selection criteria applied to the transverse momentum p^ of the spin-2 resonance are also shown. For 
the UC scenario no selection is applied. 


4 Data and simulated samples 

The data presented in this paper were recorded by the ATEAS detector during the 2012 EHC run with 
proton-proton collisions at a centre-of-mass energy of 8 TeV, and correspond to an integrated luminosity 
of 20.3 fb“^. For the H ^ yy and H ZZ* M channels, the data collected in 2011 at a centre- 
of-mass energy of 7 TeV corresponding to an integrated luminosity of 4.5 fb“\ are also used. Data 
quality requirements are applied to reject events recorded when the relevant detector components were 
not operating correctly. More than 90% of the recorded luminosity is used in these studies. The trigger 
requirements used to collect the data analysed in this paper are the same as those described in previous 
publications [18-20]. They are only briefly recalled in fhe following sections. 

The Monfe Carlo (MC) samples for fhe backgrounds and for fhe SM Higgs boson signal are fhe same as 
those used for the analyses described in Refs. [18-20], whereas new non-SM signal samples have been 
simulated. An overview of the signal samples is given in Section 4.1. 

The effects of the underlying event and of additional minimum-bias interactions occurring in the same or 
neighbouring bunch crossings, referred to as pile-up in the following, are modelled with Pythia 8 [23]. 
The ATEAS detector response is simulated [24] using either Geant 4 [25] alone or combined with a 
parameterised Geant 4-based calorimeter simulation [26]. 
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4.1 SM Higgs boson and BSM signal samples 


The SM Higgs boson ggF production for all analyses is modelled using the Powheg-Box [27] generator 
at next-to-leading order (NLO), interfaced to Pythia 8 for parton showering and hadronisation and to 
simulate multi-parton interactions. To improve the modelling of the SM Higgs boson pT. a reweighting 
procedure is applied. This procedure applies a weight depending on the px of the Higgs boson to each 
event. The weights are chosen in order to reproduce the prediction of the next-to-next-to-leading-order 
(NNLO) and next-to-next-to-leading-logarithms (NNLL) dynamic-scale calculation given by the hres2.1 
program [28, 29]. 

For the // —> yy analysis, the signal samples are generated at several values of the Higgs boson mass 
niii around 125 GeV. The samples are used to obtain a parameterisation of the signal yields and of the 
invariant mass distribution of the two-photon system as continuous functions of mn (both inclusively and 
for each category in the analysis, as described in Section 5.2). The spin-2 samples are generated using 
the MADGRAPH5_aMC@NLO [16] program with LO accuracy for zero, one, and two additional partons, 
and with subsequent matching of the matrix-element calculation with a model of the parton shower, 
underlying event and hadronisation, using Pythia 6 [30]. 

In the H ZZ* analysis the signal samples representing the production and decay of Higgs 

bosons with spin-0 and different parities are generated as follows. The SM Higgs boson production via 
gluon fusion at the mass mn - 125.5 GeV is simulated using the Powheg-Box generator. For the non- 
SM signals, the decays of the generated Higgs bosons are simulated, according to the Higgs boson parity 
assumptions, using the JHU [9, 10] MC generator at leading order (LO). The spin-2 samples are generated 
using the MADGRAPH5_aMC@NLO MC generator, as for the H ^ yy analysis. 

For the H —> WW* —> evjiv analysis, the SM Higgs boson signal is generated at mu = 125 GeV using 
the Powheg-Box Monte Carlo generator. The spin-0 BSM signal samples are generated using Mad- 
GRAPH5_aMC@NLO. The signal samples representing the production and decay of Higgs bosons with 
spin-2 are generated using the MADGRAPH5_aMC@NLO MC generator, as for the H ^ yy analysis. 

For studies of the tensor structure of the HVV decay, all simulated signal samples are obtained by using 
the matrix element (ME) reweighting method applied, as explained in the following, to a sample generated 
with non-zero values of the BSM couplings. The reweighting procedure is validated against samples 
produced at different values of the couplings, to ensure that the distributions of the CP-sensitive final- 
state observables and of their correlations are reproduced correctly. For the H ZZ* —> analysis, 
the MC production is only performed for one set of tensor couplings: g\ = f, g 2 = f + i, qa = f + i. 
All other configurations of couplings are obtained by reweighting this sample at generator level. The 
ratios of the corresponding squares of ME values calculated at EO are used as weights. To calculate these 
ME values, the JHUGenME [10] program is used. In the H WW* —> evjuv analysis, only one MC 
sample is generated, using MADGRAPH5_aMC@NEO with parameters /csm - L kaww = 2, khww - 2, 
cos(q:) = 0.3, and all other samples are obtained from it by reweighting the events on the basis of the ME 
amplitudes. 

In all the analyses presented in this paper, the mass of the Higgs boson is fixed fo 125.4 GeV [12]. 


4.2 Background samples 


The MC simulated samples for the backgrounds, as well as for the determinations of the corresponding 
cross sections, are the same as those adopted in Refs. [18-20]. In the H ^ yy analysis, the background 
is dominated by prompt yy events, with smaller contributions from y-jet events. For the H ZZ* 
analysis, the major background is the non-resonant ZZ* process, with minor contributions from the tt and 
Z-i-jets processes. For the H WW* —> ev/iv analysis, the dominant backgrounds are non-resonant W 
boson pair {WW) production, ti and single-top-quark production, and the Z/y* process followed by the 
decay to tt final states. 


5 Tests of fixed spin and parity hypotheses 

The H ^ yy and H ZZ* —> At analyses are improved with respect to the previous ATLAS publication 
of Ref. [4]. These analyses are described in some detail in the following subsections. The spin and 
parity analysis in the H —> WW* —> evjuv channel has also been improved, as discussed in detail in a 
separate publication [8]. In the following, only a brief overview of this analysis is given. The expected 
and observed results of the individual channels and of their combination are presented in Section 5.5. 


5.1 Statistical treatment 

The analyses rely on discriminant observables chosen to be sensitive to the spin and parity of the signal. 

A likelihood function, Xfdata | ,n,§), that depends on the spin-parity assumption of the signal is con¬ 

structed as a product of conditional probabilities over binned distributions of the discriminant observables 
in each channel: 


•^chann. •M)ins 

Xfdata I f, yi,e)= Y] I ^ ' ^0) - (4) 

j i 

where /ly represents the parameter associated with the signal rate normalised to the SM prediction in each 
channel j? The symbol ^ represents all nuisance parameters. The likelihood function is a product of 
Poisson distributions P corresponding to the observation of A, y events in each bin i of the discrimin¬ 
ant observables, given the expectations for the signal, \6), and for the background, R,- i{9). Some 
of the nuisance parameters are constrained by auxiliary measurements. Corresponding constraints are 
represented by the functions Jlj{0). 

While the couplings are predicted for the SM Higgs boson, they are not known a priori for the alternative 
hypotheses, defined as 7^^, as discussed in Section 3. In order to be insensitive to assumptions on the 
couplings of the non-SM resonance (the alternative hypotheses) to SM particles, the numbers of signal 
events in each channel, for each different LHC centre-of-mass energy and for each tested hypothesis, are 
treated as independent parameters in the likelihood and fitted to the data when deriving results on the spin 
and parity hypotheses. 


^ Here channel can be used to indicate different categories in the same final state when producing results for individual decay 
channels, or different final states when combining them. 
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The test statistic q used to distinguish between the two spin-parity hypotheses is based on a ratio of 
profiled likelihoods [31, 32]: 


q - log 




(5) 


where £,{J^,jujp,9jp) is the maximum-likelihood estimator, evaluated under either the SM = O'*' or 

the alternative 7^^ spin-parity hypothesis. The parameters fijp and 6jp represent the values of the signal 
strength and nuisance parameters fitted to the data under each spin and parity hypothesis. The distributions 
of the test statistic for both hypotheses are obtained using ensemble tests of MC pseudo-experiments. For 
each hypothesis test, about 70000 pseudo-experiments were generated. The generation of the pseudo¬ 
experiments uses the numbers of signal and background events in each channel obtained from maximum- 
likelihood fits to data. In the fits of each pseudo-experiment, these and all other nuisance parameters are 
profiled, i.e. fiffed fo fhe value fhaf maximises fhe likelihood for each value of fhe paramefer of inferesf. 
When generating fhe disfribufions of fhe fesf sfafislic for a given spin-parily hypofhesis, fhe expecfafion 
values of fhe signal sfrengfhs are fixed fo fhose obfained in fhe fif fo fhe dafa under fhe same spin-parify 
assumption. The disfribufions of q are used fo defermine fhe corresponding p-values pCTgjyj) = and 
For a fesfed hypofhesis 7^^, fhe observed (expecfed) p-values are obfained by infegrafing 
fhe corresponding disfribufions of fhe fesf sfafisfic above fhe observed value of q (above fhe median of fhe 
q disfribufion). When fhe measured dafa are in agreemenf wifh fhe fesfed hypofhesis, fhe observed 
value of q is disfribufed such fhaf all p-values are equally probable. 

Very small values of fhe infegral of fhe disfribufion of fhe fesf sfafisfic for fhe 7^^ hypofhesis, correspond¬ 
ing fo large values of q, are inferprefed as fhe dafa being in disagreemenf wifh fhe fesfed hypofhesis in 
favour of fhe SM hypofhesis. 

The exclusion of fhe alfemafive 7^^ hypofhesis in favour of fhe SM hypofhesis is evaluafed in ferms 
of fhe modified confidence level CLs(7^j), defined as [33]: 


CLs(7,^^) - 




1 - K^sm) 


( 6 ) 


5.2 Spin analysis in the H yy channel 

The analysis in fhe H ^ yy channel is sensifive fo a possible spin-2 sfafe. Since fhe spin-2 models 
invesfigafed in fhe presenf paper are differenl from fhose assumed in Ref. [4], fhe analysis has been 
redesigned, fo improve ifs sensifivify fo fhe new models. 

The selecfion of H ^ yy candidafe evenfs is based on fhe procedure of ofher recenf ATLAS H ^ yy 
analyses (see for example Ref. [20]). Evenfs are selecfed if fhey satisfy a diphofon frigger criferion 
requiring loose phofon idenfificafion, wifh fransverse momenfum pT fhresholds of 35 GeV and 25 GeV 
for fhe phofon wifh fhe highesf (yi) and second-highesf ( 72 ) px, respectively. During fhe offline selecfion 
fwo phofons are furfher required fo be in a fiducial pseudorapidify region, defined by \q^\ < 2.37, where 
fhe barrel/end-cap fransifion region 1.37 < \q^\ < 1.56 is excluded. The fransverse momenfum of fhe 
phofons musf safisfy p^' > 0.35 • niyy and > 0.25 • niyy, and only evenfs wifh a diphofon invarianf 
mass niyy befween 105 GeV and 160 GeV are refained. For fhe evenfs passing fhis selecfion, a furfher 
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requirement is applied on the diphoton transverse momentum, < 300 GeV, motivated by the assumed 
validity limit of the spin-2 EFT model, as explained in Section 3. After this selection, 17 220 events are 
left at a centre-of-mass energy ^[s -1 TeV and 94 540 events ai ^fs TeV. 

Kinematic variables sensitive to the spin of the resonance are the diphoton transverse momentum and 
the production angle of the two photons, measured in the Collins-Soper frame [34]: 

+ {p!^ Imyy'f 

where is the separation in pseudorapidity of the two photons. 




Figure 1: Expected distributions of kinematic variables sensitive to the spin of the resonance considered in the 
H ^ yy analysis, (a) transverse momentum of the yy system pl^ and (b) the production angle of the two photons 
in the Collins-Soper frame | cos 0*|, for a SM Higgs boson and for spin-2 particles with three different choices of 
the QCD couplings. 

The predicted distributions of these variables, for events passing the selection, are shown in Figure 1, for 
a SM Higgs boson and for a spin-2 particle with different QCD couplings. For the Kq + Kg cases, the 
enhanced high-p^^ tail offers the best discrimination, whereas for Kq - Kg the most sensitive variable is 
|cos6»*|. 

To exploit the signal distribution in both pJ^ and | cos 0*|, the selected events are divided into 11 mutually 
exclusive categories: 10 categories (labelled from Cl to CIO) collect events with <125 GeV, divided 
into 10 bins of equal size in | cos 6*\, while the 11'^ category (labelled Cll) groups all events with pl^ > 
125 GeV. As described in Section 3, for the non-UC spin-2 models the analysis is performed with two 
pY selections, namely pl^ < 300 GeV and pY <125 GeV: the latter case corresponds to not using the 
1T*’ category. 

The number of signal events above the continuum background can be estimated through a fit to the 
observed niyy distribution in each category. The niyy distribution is modelled in each category as the sum 
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of one-dimensional probability density functions (pdf) for signal and background distributions: 




[c] 
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'^^^{niyy) + + n, 


[C] 


)fs^{myy) 


bias''-'S 
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-I- n , + n, 


bias 


( 8 ) 


where J is the spin hypothesis, and are the background and the signal yield in category c, and 
fg\myy), fg‘^\myy) aTC thc Myy pdfs for the background and the signal, respectively. The signal pdf 
f^‘^\myy) is modelled as a weighted sum of a Crystal Ball function, describing the core and the lower 
mass tail, and of a Gaussian component that improves the description of the tail for higher mass values. 
For each category, f^‘^\myy) is fitted to the simulated myy distribution of the SM Higgs boson and verified 
fo be consisfenf also wifh fhe spin-2 models. The background pdf fg\myy) is empirically modelled as an 
exponential of a firsl- or second-degree polynomial. The choice of such a parameferisafion can induce a 
bias (“spurious signal”) in fhe filled signal yield, which is accounted for by fhe term The size of fhe 
expecled bias is determined as described in Refs. [20, 22], and ranges belween 0.6 and 4 evenls, depending 
on fhe category (wifh fhe signal ranging from 15 lo more lhan 100 evenls). In fhe sfalislical analysis, 
is conslrained for each calegory by multiplying fhe likelihood function by a Gaussian function cenfred al 
zero and wifh a widlh determined by fhe size of fhe expected bias. 


Defining ns as fhe folal signal yield (summed over all calegories), fhe expected fraction of signal evenls 

[c] 

n ^ 

belonging fo each category, depends on fhe spin hypofhesis J. The values of exlracled 

ns 

from fhe dala can be compared lo Iheir expected values for each spin hypofhesis, as shown in Figure 2 for 
fhe dala collected al = 8 TeV. 



Category 


2 

lo 

c 

O) 

iy) 



Cl C2 C3 C4 C5 C6 C7 C8 C9 CIO 

Category 


(a) 


(b) 


Figure 2: Observed signal fraction per category for the H ^ yy analysis, and comparison to expected values for a 
SM Higgs boson and for a spin-2 particle with different choices of QCD couplings, (a) the 11 categories described 
in the text are displayed, corresponding to the < 300 GeV selection; (b) the high-pl^ category is discarded 
and the signal fractions are renormalised over the 10 remaining categories, corresponding to the <125 GeV 
selection. 


For the non-UC scenario the iF*’ (high-p!^^) category provides strong discrimination power against the 
non-SM hypothesis, as visible in Figure 2(a). 

To discriminate between the SM spin-0 (/^j^ - 0''“) and alternative spin-2 hypotheses (T^j), two likelihood 
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functions £.jp , £.jp are built, following the general approach described in Eq. (4): 


In 


+ "bias) - Z 

c ee[c] 




where Xc runs over all categories and Yjee[c^ runs over all events in category c. The total signal yield ns is 
a free parameter in the likelihood model. The spin hypothesis being tested enters the likelihood function 
through the fractions of signal per category, . 

Several systematic uncertainties enter this model. They are implemented for each spin hypothesis as 
nuisance parameters, 9j, constrained by multiplicative Gaussian terms in the likelihood function (not 
included in Eq. (9) for simplicity). 

The signal fractions, for the SM Higgs boson are affected by uncertainties on the pj spectrum of 
the resonance and on the size of the interference between the resonance and continuum production. The 
former is computed as described in Ref. [20]. The relative impact on the signal fractions is less than ±1% 
for categories 1 to 8 < 125 GeV and |cos0*| < 0.8), and becomes as large as ±13% for categories 

10 and 11. The correction for the interference is evaluated according to Refs. [35, 36]. The systematic 
uncertainty is conservatively assumed to equal the correction itself, and its relative impact ranges between 
±0.1% and ±1.8%. 

No systematic uncertainty is assigned to the simulated p^ distribution of the spin-2 models. The effect of 
the interference between the resonance and continuum production is essentially not known, as it depends 
on the width, E^, of the resonance, which is unknown. The results presented here only hold under the 
assumption of a narrow width for the resonance, such that interference effects can be neglected. 

Additional systematic uncertainties come from the calibration of the photon energy scale and energy 
resolution and affect the signal parameterisation f^‘^\ These uncertainties are evaluated as described in 
Ref. [12]. 


5.3 Spin and parity analysis in the H -> WW* -> evfxv channel 

The analysis of the spin and parity in the H WW* —> evpv channel is described in detail in a separate 
publication [8]. In the following a brief summary is provided. The selection is restricted to events contain¬ 
ing two charged leptons of different flavour (one electron and one muon). The evpv channel is the most 
sensitive one [19]. The same-flavour channels {evev and pvpv) are not expected to add much in terms 
of sensitivity due to the presence of large backgrounds that cannot be removed without greatly reducing 
the acceptance of the alternative models considered in this analysis. The leading lepton is required to 
have pt > 22 GeV and to match the object reconstructed by the trigger, while the sub-leading lepton 
needs to have pt > 15 GeV. While the spin-0 analyses select only events with no jets in the final state 
(no observed jets with pj > 25 GeV within \r]\ < 2.5 or with px > 30 GeV within 2.5 < Ipl < 4.5), the 
spin-2 analysis enlarges the acceptance by allowing for zero or one jet (selected according to the above 
mentioned criteria). 

The major sources of background after the dilepton selection are Z/y*±jets (Drell-Yan) events, diboson 
{WW, WZly*, ZZjy*), top-quark {tt and single top) production, and W bosons produced in association 
with hadronic jets (W±jets), where a jet is misidentified as a lepton. The contribution from misidentified 
leptons is significantly reduced by the requirement of two high-px isolated leptons. Drell-Yan events 
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are suppressed through requirements on some of the dilepton variables^ > 20 GeV, < 2.8), 
while a cut on mu {mu < 80 GeV) targets the WW background. For alternative spin models with non- 
universal couplings, as discussed in Section 3, an additional upper bound is imposed on the Higgs boson 
Pt, reconstructed as the transverse component of the vector sum of the momenta of the two charged 
leptons and the missing transverse momentum. Additionally, for events containing one jet, which include 
substantial top-quark and IF-i-jets backgrounds, Z^-Jet and Z —> vetoes are applied, together with 
transverse mass requirements: the larger of the transverse masses of the two W bosons (each computed 
using the corresponding lepton and the missing transverse momentum) in the event is required to be 
larger than 50 GeV, while the total transverse mass of the WW system (defined with the two leptons and 
the missing transverse momentum) is required to be below 150 GeV. 

Control regions (CRs) are defined for fhe WW, top-quark and Drell-Yan backgrounds, which are the 
most important ones after the topological selection described above. The CRs are used to normalise 
the background event yields with a fit to the rates observed in data. The simulation is then used to 
transfer these normalisations to the signal region (SR). The IT-i-jets background is estimated entirely from 
data, while non-lTlT diboson backgrounds are estimated using MC simulation and cross-checked in a 
validation region. 

After the signal region selection, 4730 and 1569 candidate events are found in data in the 0-jet and 1-jet 
categories, respectively. For the latter category, the number decreases to 1567 and 1511 events when 
applying a selection on the Higgs boson px of less than 300 GeV and less than 125 GeV, respectively. In 
total 218 (77) events are expected from a SM Higgs boson signal in the 0-jet (1-jet) category, while about 
4390 (1413) events are expected for the total background. 

A BDT algorithm is used in both the fixed spin hypothesis tests and the tensor structure analyses. For 
spin-2 studies, the strategy follows the one adopted in Ref. [4], with the main difference being that the 
1-jet channel has been added. Two BDT discriminants are trained to distinguish between the SM hy¬ 
pothesis and the background (BDTq), and the alternative spin hypothesis and the background (BDT 2 ). 
Both BDTs employ the same variables, namely mu, p^, A0« and mj, which provide the best discrim¬ 
ination between signal hypotheses and backgrounds, also in the presence of one jet in the final state. All 
background components are used in the trainings. In total, five BDT 2 trainings are performed for the 
alternative spin hypotheses (one for the spin-2 UC scenario and two for each of the two spin-2 non-UC 
hypotheses corresponding to the different Pj selections), plus one training of BDTq for the SM Higgs 
boson hypothesis. 

For the spin-0 fixed hypothesis test and HWW tensor structure studies, the first discriminant, BDTq, is 
the same as the one used for the spin-2 analysis, trained to disentangle the SM hypothesis from the back¬ 
ground. A second BDT discriminant, BDTcp, is obtained by training the SM signal versus the alternative 
signal sample (the pure CP-even or CP-odd BSM hypotheses), and then applied to all CP-mixing frac¬ 
tions. No background component is involved in this case. The variables used for the BDTcp trainings are 
mu, ^(pu, Pj and the missing transverse momentum for the CP-even analysis and mu, ^(pu, Euw and 
Apx for the CP-odd analysis. The training strategy is different from the one used in the spin-2 analysis 
because, while the spin-2 signal is very similar to the background, the spin-0 signals are all similar to 
each other, while being different from the main background components. Therefore, in the latter case, 
training the signal hypotheses against each other improves the sensitivity. The resulting BDT variable is 

Throughout this section, the following variables are used: and niu are the transverse momentum and the invariant mass 

of the two-lepton system, respectively, A(f>te is the azimuthal angular difference between the two leptons, mj is the transverse 
mass of the reconstructed Higgs boson decay system, Apj is the absolute value of the difference befween the momenta of the 
two leptons and Grw = Pj - 0.5Pj + 0.5/)™“, where is the missing transverse momentum. 
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afterwards used in binned likelihood fits to test the data for compatibility with the presence of a SM or 
BSM Higgs boson. 

Several sources of systematic uncertainty are considered, both from experimental and theoretical sources, 
and are described in detail in Ref. [8]. The correlations induced among the different background sources 
by the presence of other processes in the control regions are fully taken into account in the statistical 
procedure. The most important systematic uncertainties are found to be those related to the modelling of 
the WW background, to the estimate of the fk+jets background (originating from the data-driven method 
employed) and, for the spin-2 results in particular, to the Z ^ tt modelling. 


5.4 Spin and parity analysis in the H -> ZZ* -> 4( channel 

The reconstruction of physics objects and event selection used for the H ZZ* Ai analysis is identical 
to the one presented in Ref. [12]. The main improvement with respect to the previous ATLAS publication 
of Ref. [4] is the introduction of a BDT discriminant designed to optimise the separation between the 
signal and the most relevant background process. 

Events containing four reconstructed leptons (electrons or muons) in the hnal state are selected using 
single-lepton and dilepton triggers. The selected events are classihed according to their final state: 
4jU, lelji, Ijjle and Ae, where for the decay modes 2 e 2 // and Ijile the first pair is defined to be the 
one with the dilepton mass closest to the Z boson mass. Each muon (electron) must satisfy pj > 6 GeV 
(Pt > 7 GeV) and be measured in the pseudorapidity range Ipl < 2.7 {\r}\ < 2.47). Higgs boson can¬ 
didates are formed by selecting two same-flavour, opposite-charge lepton pairs in an event. The lepton 
with the highest px in the quadruplet must have px > 20 GeV, and the leptons with the second- and 
third-highest px must have px > 15 GeV and px >10 GeV, respectively. The lepton pair with the mass 
closest to the Z boson mass is referred to as the leading lepton pair and its invariant mass as m\ 2 . The 
requirement 50 GeV < m \2 < 106 GeV is applied. The other lepton pair is chosen from the remain¬ 
ing leptons as the pair closest in mass to the Z boson. Its mass, denoted hereafter by ^ 34 , must satisfy 
12 GeV < m 34 < 115 GeV. Eurther requirements are made on the impact parameters of the leptons 
relative to the interaction vertex and their isolation in both the tracker and calorimeter. 

The main background process affecting the selection of ^ ZZ* —> A£ events is the non-resonant 
production of ZZ* pairs. This background has the same hnal state as the signal events and hereafter is 
referred to as the irreducible background. It is estimated from simulation and normalised to the expected 
SM cross section calculated at NEO [37, 38]. The reducible sources of background come from Z-i-jets and 
tt processes, where additional leptons arise due to misidentihed jets or heavy-havour decays. The rate 
and composition of the reducible backgrounds are evaluated using data-driven techniques, separately for 
the two hnal states with sub-leading muons it + and those with sub-leading electrons ii + ee. 

Only events with an invariant mass of the four-lepton system, denoted by m 4 />, satisfying the signal region 
dehnition 115 GeV < <130 GeV are selected. The expected signal and background yields in the 

signal region and the observed events in data are reported in Table 3. 

The choice of production and decay angles used in this analysis is presented in Eigure 3, where the 
following dehnitions are used: 

• 6 \ and 62 are dehned as the angles between hnal-state leptons with negative charge and the direction 
of hight of their respective Z bosons, in the four-lepton rest frame; 
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SM Signal 

ZZ* 

tt,Z ± jefs 

Tofal expecfed 

Observed 

Vi - 7 TeV 

An 

1.02±0.10 

0.65±0.03 

0.14±0.06 

1.81±0.12 

3 

2fi2e 

0.47±0.05 

0.29±0.02 

0.53±0.12 

1.29±0.13 

1 

2e2n 

0.64±0.06 

0.45±0.02 

0.13±0.05 

1.22±0.08 

2 

Ae 

0.45±0.04 

0.26±0.02 

0.59±0.12 

1.30±0.13 

2 

Tofal 

2.58±0.25 

1.65±0.09 

1.39±0.26 

5.62±0.37 

8 

Vi - 8 TeV 

An 

5.81±0.58 

3.36±0.17 

0.97±0.18 

10.14±0.63 

13 

2n2e 

3.00±0.30 

1.59±0.10 

0.52±0.12 

5.11±0.34 

8 

2e2n 

3.72±0.37 

2.33±0.11 

0.84±0.14 

6.89±0.41 

9 

Ae 

2.91±0.29 

1.44±0.09 

0.52±0.11 

4.87±0.32 

7 

Tofal 

15.4 ±1.5 

8.72±0.47 

2.85±0.39 

27.0 ±1.6 

37 


Table 3: Expected signal, background and total yields, including their total uncertainties, and observed events in 
data, in the 115 GeV < niiK <130 GeV signal region. The number of expected signal events is given for a SM Higgs 
boson mass of 125.5 GeV. 


• (1) is the angle between the decay planes of two lepton pairs (matched to the two Z boson decays) 
expressed in the four-lepton rest frame; 

• Oi is the angle between the decay plane of the leading lepton pair and a plane defined by the Zi 
momentum (the Z boson associated with the leading lepton pair) in the four-lepton rest frame and 
the positive direction of the collision axis; 

• 6* is the production angle of the Z\ defined in fhe four-lepfon resf frame. 

The final-sfafe observables sensifive fo fhe spin and parity of a boson decaying fo ZZ* —> A£ are fhe fwo 
production angles 6* and <l)i and fhe fhree decay angles O, 9\ and 02- In the case of a spin-0 boson, 
fhe differential producfion cross secfion does nol depend on fhe producfion variables cos(0*) and Oi. If 
should be nofed fhaf, as fhe Higgs boson mass is below 2mz, the shapes of the mass distributions of the 
intermediate Z bosons, m \2 and mj,A, are sensitive to the spin and parity of the resonance. In Figure 4 the 
distributions of the final-state observables sensitive to the spin and parity of the decaying resonance are 
presented. The distributions are shown for the SM and - 0“ simulated events, as well as for 

ZZ* production and reducible backgrounds in the signal region 115 GeV < <130 GeV. The events 

observed in data are superimposed on each plot. 

Two approaches were pursued to develop the discriminants used to distinguish between different spin 
and parity hypotheses. The first uses the theoretical differential decay rate for the final-state observables 
sensitive to parity to construct a matrix-element-based likelihood ratio analysis (7^-MELA). The second 
approach is based on a BDT. 

For the 7^-MELA approach [3, 9], the probability of observing an event with given kinematics can 
be calculated. This probability is corrected for detector acceptance and analysis selection, which are 
obtained from the simulated signal MC samples. The full pdf also includes a term for incorrect pairing 
of the leptons in the 4/r and Ae channels. Eor a given pair of spin-parity hypotheses under test, the final 
discriminanf is defined as fhe ratio of fhe pdf for a given hypofhesis fo fhe sum of fhe pdfs for bofh 
hypofheses. 


16 











Figure 3: Definitions of the angular observables sensitive to the spin and parity of the resonance in the 
X^ZZ* ^ 4( decay. 


For the BDT approach, a discriminant is formed for each pair of spin-parity states to be tested, 
by training a BDT on the variables of simulated signal events which fall in the signal mass window 
115 GeV < m^i < 130 GeV. For the O'*' versus 0“ test, only the parity-sensitive observables O, 6 \, 62 , 
m \2 and 0234 are used in the BDT training. For the spin-2 test, the production angles 6 * and <l)i are also 
included. 

Both analyses are complemented with a BDT discriminant designed to separate the signal from the ZZ* 
background. These discriminants are hereafter referred to as BDTzz- For the 7^-MELA analysis, the 
BDTzz discriminant is fully equivalent to the one described in Refs. [12, 18]. For the BDT analysis the 
discriminating variables used for the background BDTzz are the invariant mass, pseudorapidity, and 
transverse momentum of the four-lepton system, and a matrix-element-based kinematic discriminant 
Kd defined in Ref. [16]. The results from both methods are obtained from likelihood fits to the two- 
dimensional distributions of the background BDTs and of the spin- and parity-sensitive discriminants. 
In this way, the small correlation between these variables are taken into account in the analyses. The 
distribution of the background discriminant BDTzz versus the 7^-MELA discriminant is presented in 
Eigure 5 for the SM = 0'*' signal, the backgrounds, and the data. The projections of this distribution on 
the 7^-MEEA and the BDTzz variables, for different signal hypotheses, the backgrounds, and the data, 
are shown in Eigure 6 . In this paper, only results based on the 7^-MEEA approach are reported. The 
BDT approach was used as a cross-check and produced compatible results. 

Two general types of systematic effects impact the analyses using fixed spin and parity hypotheses: un¬ 
certainties on discriminant shapes due to experimental effects, and uncertainties on background normal¬ 
isations from theory uncertainties and data-driven background estimates. The systematic uncertainties on 
the shape are included in the analysis by creating discriminant shapes corresponding to variations of one 
standard deviation in the associated sources of systematic uncertainty. The systematic uncertainties on 
the normalisation are included as additional nuisance parameters in the likelihood. 

The list of sources of systematic uncertainty common to all ATEAS H ZZ* analyses is presented 
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Figure 4: Distributions of some of the final-state observables sensitive to the spin and parity of the resonance in 
the H —> ZZ* ^ signal region 115 GeV < m 4 c <130 GeV for data (points with errors), backgrounds (filled 
histograms) and predictions for two spin hypotheses (SM solid line and alternatives dashed lines), (a), (b) and (c); 
invariant masses mi 2 , mj,A and decay cos 0 \, respectively; (d), (e) and (f): 4), cos 0* and (hi, respectively. 


in Ref. [18]. The relative impact of these sources on the final separation for all tested hypotheses is evalu¬ 
ated and sources affecting the final separation (given in Section 5.5) by less than ±0.5% are neglected. 

The main sources of systematic uncertainties are related to the experimental error on the Higgs boson 
mass, the modelling of the irreducible ZZ* background, the uncertainty on the integrated luminosity and 
the experimental uncertainties on the electron and muon reconstruction. The uncertainty on the Higgs 
boson mass affects the final result since it impacts the shapes of the mi 2 , ^ 34 , cos 9\ and cos 62 variables. 
For the 7^-MELA method, the uncertainty on the estimate of the fraction of 4/r and Ae candidates with an 
incorrect pairing of leptons is also considered. This uncertainty is derived by comparing the corresponding 
prediction obtained from the Powheg and JHU MC generators for the SM hypothesis. A variation of 
± 10 % of the incorrect pairing fraction is applied to all spin and parity hypotheses. 

The influence of the main systematic uncertainties on the separation between the SM - O'*' and = 0~ 
hypotheses for the 7^-MELA analysis is presented in Table 4. The total relative impact of all system¬ 
atic uncertainties on the separation between the hypotheses (expressed in terms of numbers of standard 
deviations) is estimated to be about ±3%. 
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Figure 5; The distributions of the discriminant BDTzz versus the 7^-MELA discriminant for the SM j’’ - 0^ 
Higgs boson and for the backgrounds in the H —> ZZ* —> signal region 115 GeV < m/ii <130 GeV. 


Source of the systematic uncertainty 

Relative impact 

Higgs boson mass experimental uncertainty 

±2% 

ZZ* pdf 

±0.8% 

Muon momentum scale 

±0.7% 

Zbb —> normalisation 

±0.6% 

ZZ* scale 

±0.6% 

Luminosity 

±0.6% 

ejy resolution model (sampling term) 

±0.5% 

ejy resolution model (constant term) 

±0.5% 

Z —> i£ee normalisation 

±0.5% 

Fraction of wrongly paired Ai candidates 

±0.4% 


Table 4: Relative impact of the main systematic uncertainties on the expected separation (expressed in terms of 
numbers of standard deviations) between the SM j’’ - 0^ and - 0“ hypotheses for the H ZZ* ^ Al 
7^-MELA analysis. 


5.5 Individual and combined results 

The distributions of discriminant variables in data agree with the SM predictions for all three channels, 
and exclusion ranges for alternative spin hypotheses are derived. Some examples of distributions of the 
test statistic q (defined in Section 5.1) used to derive the results are presented in Figure 7. In this figure, 
the observed value is indicated by the vertical solid line and the expected medians by the dashed lines. 
The shaded areas correspond to the integrals of the expected distributions used to compute the p-values 
for the rejection of each hypothesis. The signal strengths per decay channel and per centre-of-mass energy 
are treated as independent parameters in each fit. Their values are compatible with the SM predictions. 
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Figure 6: Distributions of the 7^-MELA and of the BDTzz discriminants in the H —> ZZ* —> signal region 
115 GeV < m 4 i <130 GeV for the data (points with errors), the backgrounds (filled histograms), and for predic¬ 
tions for several spin and parity hypotheses. The SM hypothesis is shown by the solid line while the alternative 
hypotheses are shown by the dashed lines. The signal distributions are normalised to the signal strength fitted in 
data, (a), (b) and (c): 7^-MELA discriminants for 0^ SM vs 0^, 0^ SM vs 0^ and 0^ SM vs 2+, respectively; (d), 
(e) and (f); BDTzz discriminant for 0^ SM vs 0“, 0^ SM vs 0^ and 0^ SM vs 2+, respectively. 


The results obtained from the fit to the data, expressed in terms of p-values for different tested hypotheses 
and observed CLs for the alternative hypotheses, are summarised in Tables 5 and 6. As shown in Table 5, 
the sensitivity to reject alternative hypotheses is driven by the H ZZ* ^ 41 and the H WW* —> evfiv 
channels. The H ^ yy channel has sizeable sensitivity only to spin-2 models where the Pj < 125 GeV 
selection is not applied. In all cases the data prefer the SM hypothesis to the alternative models, with 
the exception of some of the spin-2 models for the H ^ yy channel. In this case both hypotheses have 
similar observed p-values, but neither of the two is below 10% . 

As summarised in Table 6, the p-values of the combined results for the three channels show good agree¬ 
ment between the data and the SM hypothesis for all performed tests. All tested alternative hypotheses 
are rejected at a more than 99.9% confidence level (CL) in favour of the SM hypothesis. 


6 Study of CP-mixing and of the HVV interaction tensor structure 


Following the discussion in Section 3, measurements of the HVV interaction tensor couplings /csm. kavv^ 
khvv and of the mixing angle a are performed. The measurements consist of fitting the ratios of couplings 
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q 



q 


(a) 


(b) 


H ^ ZZ* ^ 4Z 

\ls = 7 TeV, 4.5 fb"' 

\ls = 8 TeV, 20.3 fb' 

H -> WW* ^ evfiv 

•li=S TeV, 20.3 fb' 

H -> yy 

Vs = 7 TeV, 4.5 fb' 

Vs = 8 TeV, 20.3 fb' 



lO" 

10 ^ 

10 ^ 

10 

1 


ATLAS 


■ Data 

■ 0*SM 

2* (K^2Kg, p.^<125 QeV) 


H ^ ZZ* ^ 4Z 

Vs = 7 TeV, 4.5 fb' 

Vs = 8 TeV, 20.3 fb' 

H -> WW* ^ ev/iv 

Vs = 8 TeV, 20.3 fb' 

H -> yy 

Vs = 7 TeV, 4.5 fb' 

Vs = 8 TeV, 20.3 fb' 



(c) (d) 

Figure 7; Examples of distributions of the test statistic q defined in Section 5.1, for the combination of decay 
channels, (a); 0^ versus 0“; (b): 0^ versus 0^; (c): 0^ versus the spin-2 model with universal couplings (k^ - Kg)-, 
(d); 0^ versus the spin-2 model with Kg - iKg and the px selection at 125 GeV. The observed values are indicated 
by the vertical solid line and the expected medians by the dashed lines. The shaded areas correspond to the integrals 
of the expected distributions used to compute the p-values for the rejection of each hypothesis. 
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H ^ yy 


Tested Hypolhesis 

alt 

^ exp,/i=l 

Pexp.ft=fi 

pTs 

PoL 

Ohs. CEs (%) 

2+iKq = Kg) 

0.13 

7.5 -10-2 

0.13 

0.34 

39 

2* (Kg = 0; pt< 300 GeV) 

4.3 • 10-4 

< 3.1 - 10-2 

0.16 

2.9-10-4 

3.5-10-2 

2* (Kg = 0; pt < 125 GeV) 

9.4 • 10-2 

5.6-10-2 

0.23 

0.20 

26 

2+(Kg = 2Kg\ Pi < 300 GeV) 

9.1 • 10-4 

< 3.1 - 10-2 

0.16 

8.6-10-4 

0.10 

2-*(Kg = 2Kg\ Pi < 125 GeV) 

0.27 

0.24 

0.20 

0.54 

68 



H ^ 

WIV* 

-> evpv 


Tested Hypolhesis 

alt 

^ exp,/i=l 



ptL 

Ohs. CEs (%) 

% 

0.31 

0.29 

0.91 

2.7-10-2 

29 

0” 

6.4-10-2 

3.2-10-2 

0.65 

1.2-10-2 

3.5 

2+(Kg ^ Kg) 

6.4-10-2 

3.3-10-2 

0.25 

0.12 

16 

2+(Kg = 0; Pi < 300 GeV) 

1.5-10-2 

4.0-10-2 

0.55 

3.0-10-2 

0.6 

2^ (Kg = 0; pt < 125 GeV) 

5.6-10-2 

2.9-10-2 

0.42 

4.4-10-2 

7.5 

2+(Kg = 2Kg\ Pi < 300 GeV) 

1.5-10-2 

4.0-10-2 

0.52 

3.0-10-2 

0.7 

2-*(Kg = 2Kg\ Pi < 125 GeV) 

4.4-10-2 

2.2-10-2 

0.69 

7.0-10-2 

2.2 



H - 

ZZ* 

^ At 


Tested Hypolhesis 

alt 

^ exp,/i=l 

Pexp,fi=fi 


PoL 

Ohs. CU (%) 

K 

3.2 - 10-2 

5.2 -10-2 

0.80 

3.6 - 10-4 

0.18 

0- 

8.0 - 10-2 

3.6 -10-4 

0.88 

1.2- 10-2 

1.0-10-2 

2'*(Kg = Kg) 

3.3 - 10-2 

5.7 -10-4 

0.91 

3.6 - 10-2 

4.0-10-2 

2-* (Kg = 0; pt < 300 GeV) 

3.9 - 10-2 

9.0 -10-2 

0.95 

2.7 - 10-2 

5.4-10-2 

2-^ (Kg = 0; Pt < 125 GeV) 

4.6 - 10-2 

1.1 -10-2 

0.93 

3.0 - 10-2 

4.3-10-2 

2-*(Kg = 2Kg\ Pi < 300 GeV) 

4.6 - 10-2 

1.1 - 10-2 

0.66 

3.3 - 10-2 

0.97 

2-*(Kg = 2Kg\ Pi < 125 GeV) 

5.0 - 10-2 

1.3 -10-2 

0.88 

3.2 - 10-4 

0.27 


Table 5: Expected and observed p-values for different spin-parity hypotheses, for each of the three channels H 
yy, H ZZ* ^ and H —» WW* ev/jv. The observed CLj for the alternative hypotheses are reported 
in the last column. The expected and observed p-values and the observed CLj are defined in Section 5.5 and the 
alternative hypotheses are those described in Section 3. 


(^avv/z^sm) • tana and khvvIksm to the discriminant observables for the H —> WW* —> ev//y and H 
ZZ* —> At processes and in their combination. In the fitting procedure only one ratio of couplings 
(i^Avy/^SM) • tan a or khvvIksm is considered at a time, while the other one is assumed to be absent. 


6.1 Statistical treatment 

The measurement of the tensor structure of the HVV interaction is based on a profiled likelihood [31, 32] 
fhaf confains fhe discriminanf observables sensifive fo fhe EFT couplings. The signal rales in Ihe differenl 
channels and for differenl cenlre-of-mass energies are freafed as independenl paramefers. Therefore, 
the global signal normalisation is not used to constrain the EFT couplings. The ratios of the BSM to 
SM couplings, khvvIksm and QcavvIksm) • tan a, are each separately fit to the discriminant observables in 
data. The test statistic used to derive the confidence intervals on fhe paramefers of inleresl is q' - -2 ln(/l), 
where A is fhe profiled likelihood [31, 32]. The resulls presented in fhe following rely on fhe asympfolic 
approximafion [31, 32] for fhe lesl slafisfic. This approximation was cross-checked wifh Monte Carlo 
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Tesfed Hypofhesis 

alt 

^ exp,/ 2 =l 

Pexp,ii=fi 


^obs 

Obs. CLs (%) 

0 ; 



2.5 • 10-^ 

4.7 •10-3 

0.85 

7.1 • 

10-3 

4.7 • 

10-^ 

0“ 



1.8 • 10-3 

1.3 • 10-"^ 

0.88 

< 3.1 

• 10-3 

< 2.6 

• 10-3 

-- 

--Kg) 


4.3 • 10-3 

2.9 • 10-^ 

0.61 

4.3- 

10-3 

1.1 • 

10-3 

2*{Kq -- 

= 0; pt < 

300 GeV) 

< 3.1 • 10-3 

< 3.1 • 10-3 

0.52 

< 3.1 

• 10-3 

< 6.5 

• 10-3 

T{Kq -- 

= 0; pt < 

125 GeV) 

3.4 • 10-3 

3.9 • 10-^ 

0.71 

4.3- 

10-3 

1.5- 

10-3 


= '^Kg-, PT 

< 300 GeV) 

< 3.1 • 10-3 

< 3.1 • 10-3 

0.28 

< 3.1 

• 10-3 

<4.3 

• 10-3 

-- 

= '^Kg-, PT 

< 125 GeV) 

7.8 • 10-3 

1.2- 10-3 

0.80 

7.3- 

10-3 

3.7- 

10-3 


Table 6 : Expected and observed p-values for different spin-parity hypotheses, for the combination of the three 
channels: ^ 77 , // —> ZZ* ^ A{ and H WW* ev/iv. The observed CLs for the alternative hypothesis is 

reported in the last column. The expected and observed p-values and the observed CL, are defined in Section 5.5. 
The definitions of alternative hypotheses are given in Section 3. 


ensemble tests that confirm its validity in the range of the parameters for which the 95% CL limits are 
derived. 


6.2 Tensor structure analyses in the H WW* evfxv channel 

The H WW* eviiv analysis used to study the spin-0 tensor structure is already described in Sec¬ 
tion 5.3 and detailed in Ref. [ 8 ]. Only the 0-jet category is considered and the BDTq and BDTcp are used 
as discriminant variables in the likelihood defined fo measure fhe spin-0 fensor strucfure couplings. The 
only difference wifh respecf fo fhe spin hypofhesis fesf is fhaf, in fhis analysis, fhe BSM spin-0 couplings 
are freafed as confinuous variables in fhe fesf sfafisfic. 


6.3 Tensor structure analyses in the H ZZ* 4{ channel 

To allow for a cross-check and validafion of fhe obfained resulfs, fwo differenl fitting mefhods based on 
fhe analyfical calculafion of fhe leading-order mafrix elemenf of fhe H ZZ* ^ Ai process are used. 

The mefhod of fhe mafrix-elemenf-observable fif is based on modelling fhe disfribufions of fhe final-sfafe 
observables in each bin of coupling ratios using Monfe Carlo simulafion. Using fhe Lagrangian defined 
in Eq. (1), which is linear in fhe coupling consfanfs A'sM^ khvv and kavv ^ the differenlial cross secfion 
af each poinf in fhe phase space can be expressed as a term corresponding fo fhe SM amplifude, plus 
fwo addifional ferms, linear and quadrafic in fhe coupling consfanfs. In fhis way if is possible fo define 
fwo observables for each coupling, fhe so-called firsf- and second-order opfimal observables, upon which 
fhe amplifude depends af each poinf of fhe phase space. For each evenf, fhey confain fhe full kinemafic 
informafion abouf fhe couplings, which can fhus be exfracfed from a fif fo fheir shapes. More defails of 
fhe mefhod can be found in Refs. [39-42]. 

The observables sensitive fo fhe presence and sfrucfure of khvv and kavv considered in fhe currenf 
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analysis are defined as follows: 


Oi{khvv) = 
OliKHVv) = 
0\{KAvv,a) = 
OiinAVV, a) = 


2^[ME(^SM?^Qi khvv,i^avv=0''> Q'=Q)*-ME(/(://vv?^0; a'sm^.^avv'— Qi Q^=0)] 

|ME(yfsM?^0; khvv,/<avv=0-, o'=0)p ’ 

\ME{khvv*0\ ksu,kavv=0\ Q'=0)P 

|ME(*:sm^0; khvv,kavv=0', Q'=0)P ’ 

2^[ME(/csm^Q^ f<HVVif<Avv=0\ Q'=Q)*-ME(a:^v/\/:;^ 0; Qf=;r/2)] 

|ME(yfsM?^0; khvv,/<avv=0; o'=0)p 
ksm,i^hvv= 0', a=n/2)f 

|ME(fcsM^0; khvv,kavv=^', Q'=0)P 


( 10 ) 


Here ME(/csm» k^hvv^ kavv^ o') denotes the leading-order matrix element of the H ZZ* process. 
These definitions correspond to the first- and second-order optimal observables for a BSM amplitude with 
a three-component structure. 

The observables Oi^ii^HVv) and 0\^2{kavv^ o') are used for the khwI^sm and (^^yy/zcsM)'tan a individual 
fits respectively. In order to suppress the ZZ* background, a kinematic BDT discriminant similar to those 
described in Section 5.4 is used as an additional observable in all fits. The BDT training is performed 
independently for each final state using observables with small sensitivity to parity: 774 ^, / 7 T, 4 r> cos( 0 *) 

and Oi. This BDT discriminant is denoted hereafter by BDT(ZZ). 

To simplify their use in the analysis, all observables defined in Eq. (10) undergo a pdf transformation 
such that each observable becomes normally distributed in the Standard Model case. These transformed 
observables are referred to hereafter as TOi^ixHVv) and T0\2{kavv, ct) respectively. The distributions of 
transformed observables for the Monte Carlo signal samples generated with {khvvIk$m = 0, ±\-,kavv - 
0) and ((kAvy/zcsM) • tan a = 0, ±5; khvv = 0) are shown in Eigure 8 . The contributions of all backgrounds 
considered in this analysis are also included. By construction the TO 2 observables are sensitive to the 
modulus of the khvvI and (kavv/^sm) ■ tan a ratios: their distributions change with the strength of the 
respective coupling. These observables are insensitive to the relative sign of khvv and kavv with respect 
to a:sm- The sign sensitivity comes from the TOi observables, which are based on the interference terms: 
their distributions feature pronounced sign-dependent asymmetries. It was also found that the observables 
TO\{khvv) and T02(khvv) are linearly correlated. To maximise the population of analysis histograms 
with currently available Monte Carlo event samples, it is desirable to reduce this correlation. This is 
achieved by considering the modified observables TOi{khvv) + T02{khvv) and T0\{khvv) - T02{khvv) 
in the current analysis. 

The analysis is performed in several steps. Eirst, multi-dimensional histograms of observables are created 
in 81 bins of khvvIk^sm and (kavv/ksm) ■ tana for all fits. The predicted shapes of the observables for the 
signal are produced by reweighting the base Monte Carlo sample described in Section 4. The correspond¬ 
ing weights are derived using the analytical calculation of the H ZZ* —> 4€ matrix elements at leading 
order in perturbative QCD. The weights are calculated and applied at the Monte Carlo generator level. 

The observables used in the analysis are evaluated after detector simulation, accounting for the detector 
acceptance, resolution and reconstruction efficiency. The distributions of observables for backgrounds 
are estimated using Monte Carlo (for the irreducible background) and data-driven techniques (for the 
reducible backgrounds) described in Section 5 and Refs. [12, 18]. 

The distributions of observables are three-dimensional: T 0\{kavv^ »)> T 02{kavv, «)> BDT(ZZ) and T 0\{khvv) 
+ T02{khvv)^ T0\{khvv) - T02{khvv)^ BDT(ZZ) respectively. To obtain a reliable description for bins 
with an insufficient number of Monte Carlo events, the Kernel Density Estimation [43] smoothing pro¬ 
cedure is applied to signal and background multi-dimensional histograms. In the smoothing procedure 
the smearing is done separately in four bins of BDT(ZZ), preserving the original normalisation. 
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rfs = 8 TeV, 20.3 fb'' 


ATLAS 

H-^ZZ*-^4I 


= 0 , a = 0 


■ ^ ’ ^AW ■ 
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SM: Kg„ = 1, k,,„ = 0 , ic.yy = 0 , a = 0 


(e) 

Figure 8; Distributions of the observables used in the matrix-element-observable fit. (a); T0\{khvv) + T02{khvv), 
(b): TOi{khvv) - T02(khvv), for the Monte Carlo signal generated with (khvvIksm - 0,+1;a-avv = 0). (c): 
TOi{kavv, a), (d); T02(i<Avv, a) for the Monte Carlo signal generated with (fevv/zcsM) ■ tan a = 0, +5; khvv - 0). 
(e); BDT(ZZ) for the Monte Carlo signal generated with (khvvIksm - 0, ±\\kavv - 0). The expected background 
contributions are shown as filled histograms on each plot. 
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The final pdfs used in the fits are obtained by applying linear histogram interpolation between the multi¬ 
dimensional bins of khvvIksm and {kavvIksm) • tana. The individual likelihood functions per centre-of- 
mass energy ( ^fs) and final state (FS) are: 



KhVV Kavv 
/CSM ’ /<^SM 


tan a, 9 


n 


a, 


^KHVV KAVV ^ 7)1 

Si I-,-tan a,9\ + bi{9) 

/CSM /<^SM 


(11) 


where P is the probability density function for the data vector Q, given the signal model s and back¬ 
ground model b. The index i runs over all the bins of multi-dimensional histograms of observables and 
9 represents the vector of nuisance parameters corresponding to systematic uncertainties. Fits to data are 
performed by minimising the negative log-likelihood function with respect to the ratios of the couplings: 


L 



KhVV KAVV 
/CSM ’ /CSM 


tan a, 9 


-2innn^r 

VJ FS 1 


T<HVV KAVV 
/<^SM /<^SM 


tan a, 9 


( 12 ) 


The test statistic q' - -2 ln(/l) is defined as the profiled value of L of Eq. (12). To ensure the correctness 
of the statistical treatment and the absence of significant biases, a series of tests were performed before 
applying the fit to the data. Asimov datasets [31, 32] created from independently generated Monte Carlo 
samples with khvvIksm and {kavvIksm) • tana equal to 0, ±2, ±4, ±6, ±8 and ±10 were injected into the 
analysis procedure. The tests were repeated for samples corresponding to 1 and 100 times the LHC Run-I 
integrated luminosity. In all cases the fitted values of coupling constants were found to be in agreement 
with the injected values within statistical uncertainties. 

The results of the matrix-element-observable fit were validated and cross-checked using a nine-dimensional 
matrix-element method (9D fit). The method implements a multivariate per-event extended likelihood that 
is sensitive to both the khvvIksm and {kavvIksm) • tana mixing parameters and is based on nine exper¬ 
imental observables. The probability model is constructed with separate components for signal, the SM 
ZZ* background and the reducible background. The background components are assumed to be inde¬ 
pendent of the Higgs boson tensor structure, so all of the sensitivity to mixing parameters comes from the 
signal component. Each component depends on nine experimental observables: Pt^ 4 £, qnt, cos 9*, 

cos 9\, cos 92, O, mi 2 and (described in Section 5.4). 

The main sources of systematic uncertainty for the tensor stmcture measurements are the same as dis¬ 
cussed in Section 5 since they are based on the same four-lepton variables. Several additional sources 
of uncertainty, specific to each of the methods, are also taken into account. Eor the matrix-element- 
observable fit, the uncertainty related to the Kernel Density Estimation smoothing procedure applied to 
signal and background multi-dimensional histograms is considered. To estimate the influence of this un¬ 
certainty on the final result, a procedure similar to the one described in Section 5 is employed. The impact 
of the different sources of systematic uncertainty on the final results is evaluated by comparing the BSM 
exclusion limits obtained with a specific systematic uncertainty included or excluded in the fit, while 
excluding all other systematic uncertainties. A similar conclusion holds in the fixed hypothesis test: the 
systematic uncertainties have a very limited impact on the final result. The most important uncertainties 
are related to the estimates of the reducible backgrounds. The relative impact of these uncertainties on the 
final 95% CL exclusion limit on BSM couplings was found to be around ±1%. The second most import¬ 
ant group of sources of systematic uncertainty is related to the theoretical uncertainties on the production 
cross section of the ZZ* background process. Their relative impact on the final result is found to be less 
than ±1%. The precision of the tensor structure analysis is thus dominated by the statistical errors. 
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In this paper, only results based on the matrix-element-observable approach are reported. The 9D ap¬ 
proach was used as a cross-check and produced results compatible with the matrix-element approach. 


6.4 Individual and combined results 

The results of the tensor structure analyses performed in the H —> WW* —> evfiv channel are reported in 
Ref. [8] and, for completeness, they are also summarised in Table 7. 


Coupling ratio 

Best-fit value 

95% CL Exclusion Regions 

H ITIT* ^ evjiv 

Observed 

Expected Observed 

KnvvlKsyi 

-1.3 

[-1.2, -0.7] (-CX,, -2.2] U[-l, -0.85] U[0.4, cx,) 

{kavvIksm) ■ tana 

-0.2 

n.a. (-O 0 ,-6] 1J[5, oo) 


Table 7; Fitted values of khvvIksm and (kavvIksm) ■ tana and 95% CL excluded regions obtained in // —> WW* —> 
ev/uv analysis. The expected values are estimated for the signal strength measured in data and assuming best-fit 
values for all other nuisance parameters. Only data collected at -yT = 8 TeV are used. The symbol "n.a." denotes 
the absence of 95% CL sensitivity. 

The distributions of the test statistic for fits of khwI^sm and {kavvIksm) • tana measured in the H —> 
ZZ* —> Al analysis are shown in Figure 9. The expected curves are calculated assuming the SM 




(a) (b) 

Figure 9: Expected and observed distributions of the test statistic for hts of (a) khvvIk^u and (b) {kavvIksm) ■ tan a 
for the H ZZ* —> A{ analysis. The expected curves are calculated assuming the SM j’’ - 0^ signal and produced 
with the SM signal strength n - \ and with the signal strengths fitted to data. The horizontal dotted black lines 
represent the levels of -2 In A above which the values of coupling ratios under study are excluded above 68% and 
95% CL, respectively. 

signal, both with the SM signal strength, fi = I, and with the signal strength fitted to data, fi. The 
fitted values of khwI^sm and {kavvIksm) • tana, together with the intervals where these couplings are 
excluded at above the 95% CL, are reported in Table 8. The fitted values agree with the SM predictions 
within uncertainties. 
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Coupling ratio 

Best-fit value 

95% CL Exclusion Regions 

H ^ZZ* ^ M 

Observed 

Expected 

Observed 

khvvIksm 

-0.2 

(-00,-0.75] U[6.95,00) 

(-00,-0.75] U[2.45,00) 

{kavvIksm) ■ tana 

-0.8 

(-CX,,-2.95] U[2.95, cx,) 

(- 00 ,-2.85] U[0.95, oo) 


Table 8; Expected and observed best-fit values of khvvIksm and (kavvIksm) ■ tana and 95% CL excluded regions 
obtained in the H —» ZZ* —> analysis. The expected values are estimated for the signal strength measured in 
data and assuming best-fit values for all other nuisance parameters. The data for ■\fs - 1 TeV and Vs = 8 TeV are 
combined. 


The measurements from the H WW* —> ev/tv and H —> ZZ* channels are combined under the 
assumption that the BSM ratios of couplings khwIk^sm and QcavvIksm) • tana are the same for the W 
and Z vector bosons. A common test statistic is obtained by combining the profiled likelihoods of the 
individual channels. The expected distributions of the likelihoods, for the signal strength values obtained 
from the fits to the data (jj. = p), are presented in Figure 10. The observed distributions of profiled 


_c 




(a) (b) 

Figure 10: Expected distributions of the test statistic for the combination of H ^ WW* evfiv and H —> ZZ* ^ 4^ 
analyses as a function of BSM coupling ratios (a) khvvIksm and (b) (kavv / ksm) ■ tana. The expected values are 
estimated for the signal strengths measured in data and assuming best-fit values for all other nuisance parameters. 
The 68% and 95% CL exclusion regions are indicated as lying above the corresponding horizontal lines. The 
individual distributions for H —> WW* evfiv and H —> ZZ* ^ channels are shown. 

likelihoods for the combination ofH^ WW* —> cv/rv and H ZZ* —> measurements are presented 
in Figure 11. The asymmetric shape of the expected and observed limits in the khvvIk$m results is mainly 
due to the interference between the BSM and the SM contributions that gives maximum deviation from 
the SM predictions for negative relative values of the BSM couplings. Here the signal normalisations are 
treated as independent nuisance parameters of the different decay channels and the different centre-of- 
mass energies. The other nuisance parameters related to the experimental and theoretical uncertainties 
are treated as correlated when appropriate. The resulting 95% CL exclusion regions for the combinations 
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(a) 


(b) 




(c) (d) 

Figure 11: Expected and observed distributions of the test statistic for H —> WW* —> ev/rv and H —> ZZ* ^ 
analyses and their combinations. The distributions are shown as a function of the BSM coupling ratios khvvIk^u 
and {kavvIk^u) ■ tana. The 68% and 95% CL exclusion regions are indicated as lying above the corresponding 
horizontal lines, (a) and (b): individual H —» WW* —» evfiv, H ZZ* —» and combined observed distributions, 

(c) and (d): expected and observed combined distributions. The expected distributions are presented for the SM 
signal strength ju - I and for the signal strengths obtained from the ht to data. 


29 





















































of H ^ WW* evfiv and H ZZ* 4£ channels are listed in Table. 9. 


Coupling ratio 

Best-fit value 

95% CL Exclusion Regions 

Combined 

Observed 

Expected 

Observed 

KhvvIksm 

-0.48 

(-C50,-0.55] U[4.80, oo) 

(- 00 ,-0.73] U[0.63, oo) 

{kavvIksm) ■ tana 

-0.68 

(- 00 ,-2.33] U[2.30, oo) 

(-oo,-2.18]U[0.83,oo) 


Table 9: Expected and observed best-fit values of (a) khvvIksm and (b) {kavvIksm) ■ tana and 95% CL excluded 
regions obtained in the combination of H ^ ZZ* —> At and H WW* ev/uv analyses. The expected values are 
estimated for the signal strengths measured in data and assuming best-fit values for all other nuisance parameters. 
The signal strengths are treated independently per decay channel and per collision energy. 


7 Conclusion 

Studies of the spin and parity of the observed Higgs boson in the H ZZ* M, H ^ WW* —» evjiv 
and H ^ yy decay processes are presented. The investigations are based on 4.5 fb“^ and 20.3 fb“^ 
of pp collision data collected by the ATLAS experiment at the LHC at = 7 TeV and = 8 TeV, 
respectively. The SM Higgs boson hypothesis, corresponding to the quantum numbers = O'*', is tested 
against several alternative spin and parity models. The models considered include non-SM spin-0 and 
spin-2 models with universal and non-universal couplings to quarks and gluons. The combination of 
the three decay processes allows the exclusion of all considered non-SM spin hypotheses at a more than 
99.9% CL in favour of the SM spin-0 hypothesis. 

The tensor structure of the HVV interaction in the spin-0 hypothesis is also investigated using the H 
ZZ* and H WW* —» evpv decays. Only one BSM tensor coupling is investigated at a time, while 
the other one is set to zero. The observed distributions of the variables sensitive the ratios of the BSM to 
SM tensor couplings, khvvIksm and (k^w/A^SM) • tana, are compatible with the SM predictions. 

Values of the BSM tensor couplings outside of the intervals -0.75 < khvvIksm < 2.45 and -2.85 < 
ii<AVV/ ksm) ■ tana < 0.95 are excluded at the 95% CL for the H ZZ* 4€ process. For the H 
WW* evyuv process the ranges -2.2 < khvvIksm < -1-0 and -0.85 < khvvIksm < 0.4 and -6.0 < 
(kAyy//<'SM) • tana < 5.0 are excluded at the 95% CL. 

The results from the 7/ —> WW* evpv and H ZZ* 47 decay channels are combined under the 
assumption that the kevvIksm and (kAyy/xsM) • tana couplings have the same values for the HWW and 
HZZ processes. As a result of this combination, the regions outside of -0.73 < kevvIksm < 0.63 and 
-2.18 < {kavvIksm) • tana < 0.83 intervals are excluded at the 95% CL. The corresponding expected 
not-excluded intervals at the 95% CL, assuming the SM Higgs boson hypothesis and the signal strength 
values measured in data, are -0.55 < kevvIksm < 4.80 and -2.33 < {kavvIksm) • tana < 2.30. 
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Appendix A 


To compare the exclusion limits obtained in this analysis to other existing studies, the final results of 
this analysis are also expressed in terms of effective cross-section fractions {fg 2 , (pgi) and {fg^, (f>g 4 ). The 
definitions proposed in Section 11.4.2 of Ref. [3] and Section II of Ref. [44] are used: 


fgi - 


\gi?o-i 


l^lpcn + 1021^0-2 + 1041^0-4 


, (pi = arg 


\gi 


(13) 


Here the symbols gi, g 2 and g^ denote the SM, BSM CP-even and BSM CP-odd tensor couplings of 
the HVV scattering amplitude, respectively. The numeric coefficients cri, o ‘2 and cr 4 are effective cross 
sections of the HVV interactions calculated when only the gi-, g 2 - or 94 -related terms are present in the 
amplitude, respectively, such that gi = \,gii,j = 0 . 

When, in addition to the SM term, only one CP-even or CP-odd BSM contribution is present, the con¬ 
version between the parameterisation used in this analysis and the {fgi, tpgi) parameterisation is given by 
Eq. (13) rewritten in the following way: 


/ft “ 






(/ = 2,4), 


(14) 


where r 4 i and r 2 \ are chosen such that: 


'21 


(^HVV I ^HW 
O'SM \ ^SM 


and 


Caw 

C^SM 



tan^ a. 


(15) 


The numeric coefficients ctsm, cthvv and ctavv are effective cross sections of the HVV interaction calcu¬ 
lated when only each of the atsm-. khvv- and /fAVV-related terms is present in the Lagrangian. 

For consistency with previous measurements reported in Ref. [5], the expected and observed results of 
the current analysis of the H —> WW* evfiv and H ZZ* ^ Ai channels and for their combin¬ 
ation are expressed in terms of fgi and fgi parameters for the H ZZ* decay, (/^^^0^f) and 

{fg^, 0^4 )• These parameters are denoted hereafter by {fg 2 , (pg 2 ) and (/^ 4 , (pg 4 ). The corresponding results 
are presented in Tables 10 and 11. To obtain these results, the effective cross sections ctsm, cthvv and 
(Taw of the HZZ interaction are calculated using the MADGRAPH5_aMC@NLO Monte Carlo generator 
[16] at leading order. The ratios of cross sections used in the calculation are: cthvvIcsm = 0.349 and 
cavv/csm - 0.143, respectively. 
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Expected 95% CL limits 

H lElE* ^ ev^v 

n.a. for 0^2 - 0 

and 

fg2 <0.15; fg2 > 0.33 for (pg 2 = n 

n.a. for 0^4 - 0 

and 

n.a. for (pg 4 = n 



ZZ* Ai 

fg 2 < 0.94 for 0^2 ^ 0 

and 

fg 2 < 0.16 for (pg 2 = n 

7^4 < 0.56 for (pg 4 . = 0 

and 

fg 4 < 0.56 for 0 g 4 = TT 

Combination of H 

^ZZ* 

^ A€ and H lElE* ^ ev/iv 

fg2 < 0.89 for 4>g2 ^ 0 

and 

fg 2 < 0.096 for 0^2 ^ TT 

fg 4 < 0.43 for (pgA = 0 

and 

fgA < 0.44 for (t)g 4 = n 


Table 10: Expected limits on {fg2,4’g2) (/g4,0g4) parameters defined in Ref. [3] obtained in the analyses of 

the H —» WW* evfiv and H —> ZZ* ^ 4{ channels and for their combination. The symbol "n.a." denotes the 
absence of 95% CL sensitivity. 


Observed 95% CL limits 




H ■ 

^ ITIT* ^ evfiv 



fgl 

< 0.053 for (()g2 - 

-0 

and 

fg2 < 0.20; 0.26 < fg2 < 0.63 for <pg2 = n 

fg 4 

< 0.78 for 4)g4 - 

0 

and 

fg4 < 0.84 for 

4 >gA 

= n 

H^ZZ* ^ A£ 

fgl 

< 0.68 for (pg2 - 

0 

and 

fg2 <0.16 for 

(f>gl 

= n 

fgA 

< 0.11 for 0^4 = 

0 

and 

fg4 < 0.54 for 

4 >gA 

= n 


Combination of H ^ 

ZZ* ^ Ae and H WW’ 

Ic _^ 

evfiv 

fgl 

< 0.12 for (pg2 - 

0 

and 

fg2 <0.16 for 

(pgl 

= n 

fgA 

< 0.090 for (f>g4 - 

:0 

and 

fg4 < 0.41 for 

(pgA 

= n 


Table 11: Observed imits on (fg 2 , (pg 2 ) and (/^ 4 , (pgi) parameters defined in Ref. [3] obtained in the analyses of the 
H WW* —> ev/uv and H ZZ* —> 4{ channels and for their combination. 
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M. Campanelli^^, A. Campoverde^^^, V. A. Canepa^^^^, M. Cano Bret^^'^, J. Cantero*^, 

R. Cantrill^^^^, T. Cao^°, M.D.M. Capeans Garrido^°, I. Caprini^^®, M. Caprini^^®, M. Capua^’^’^^'’, 

R. Caputo^^, R. Cardarelli^^^^, E. Cardillo^^, T. Carli^°, G. Carlino'’^^^, E. Carminati^^^’^^'’, S. Caron^®^, 
E. Carquin^^®, G.D. Carrillo-Montoya^, J.R. Carter^^, J. Carvalho^^^®’^^®‘^, D. Casadei^*, M.P Casado^^, 
M. Casolino^^, E. Castaneda-Miranda^^^'’, A. Castelli^*^’, V. Castillo Gimenez^^^, N.E. Castro^^^®’^, 

P Catastini^’, A. Catinaccio^®, J.R. Catmore^^^, A. Cattai^°, J. Caudron^^, V. Cavaliere^^^, D. Cavalli^^®, 
M. Cavalli-Sforza^^, V. Cavasinni^^^®’^^^'’, E. Ceradini^^^®’^^^'’, B.C. Cerio^^, K. Cerny^^^, 

A.S. Cerqueira^^’’, A. Cerri^^^, E. Cerrito’^, E. Cerutti^^, M. Cerv^°, A. Cervelli^^, S.A. Cetin^^‘^, 

A. Chafaq^^^®, D. Chakraborty^®^, I. Chalupkova^^^, P. Chang^^^, J.D. Chapman^*, D.G. Charlton*^, 

C. C. Chau'^^, C.A. Chavez Barajas^^^, S. Cheatham^^^, A. Chegwidden^*^, S. Chekanov^, 

S. V. Chekulaevi59a, q.A. Chelkov^^.ft^ ^ ^ Chelstowska^^ C. Chen^\ H. Chen^^, K. Chenl4^ 

E. Chen33d.', 5 chen^^^^, X. Chen^^f, Y. Chen^^, H.C. Cheng^^ Y. Cheng^i, A. Cheplakov^^ 

E. Cheremushkina*^°, R. Cherkaoui El Moursli'^^®, V. Chernyatin^^’*, E. Cheu^, E. Chevalier^^^, 

V. Chiarella^^, J.T. Childers^ G. Chiodini^^a^ ^ S Chisholmi^ R.T. Chislett^^ A. Chitan^^^ 

M.V. Chizhov^^ K. Choi^\ S. Chouridou'^, B.K.B. Chow^°°, V. Christodoulou^^ 

D. Chromek-Burckhart^*^, J. Chudoba^^’, A.J. Chuinard^’, J.J. Chwastowski^^, E. Chytka^^^, 

G. Ciapettii32a,i32b^ ^ Ciftci^^ D. Cinca^^, V. Cindro’^ I.A. Cioara^i, A. Ciocio^^ E. Cirottol°4^’l°^^ 
Z.H. Citron^^^ M. Ciubancan^^^ A. Clark^^, B.E. Clark^^, PJ. Clark^^, R.N. Clarkei^, W. Clelandi^s, 

C. Clementl4^“’^^'^^ Y. Coadou^^ M. Coball6^“’^^4^ A. Coccaro^^^ J. Cochran^^, E. Coffey^^, 

J. G. Cogani43, B. Cole^^, S. ColelO^ A.P Colijn^O’, J. Collot^^ T. Colombo^^^ G. Compostellai°i, 

P. Conde Muinoi 26 a.i 26 b^ g Coniavitis^^ S.H. Connell^^^b^ j ^ Connelly^’, S.M. Consonni‘^^^’‘^l^ 

V. Consorti^^ S. Constantinescu^®^ C. Conta^ 2 ia,i 2 ib^ q Conti^o, E. Conventii°^^’3, M. Cooke^^ 

B. D. Cooper’^, A.M. Cooper-Sarkar^^*^, T. Cornelissen^^^, M. Corradi^*^®, E. Corriveau^^’^, 

A. Corso-Radu^^^, A. Cortes-Gonzalez^^, G. Cortiana^*^^ G. Costa^^®, M.J. Costa^^^, D. Costanzo^^^, 

D. Cote*, G. Cottin^*, G. Cowan^’, B.E. Cox*^, K. Cranmer^^°, G. Cree^^, S. Crepe-Renaudin^^, 

E. Crescioli*°, W.A. Cribbs^^^®’^^^'’, M. Crispin Ortuzar*^°, M. Cristinziani^\ V. Croft^°^, 

G. Crosetti*’^^’*’'’, T. Cuhadar Donszelmann^^^, J. Cummings^^^, M. Curatolo^^, C. Cuthbert^^**, 

H. Czirr^^\ P. Czodrowski^, S. D’Auria^^, M. D’Onofrio^^, M.J. Da Cunha Sargedas De Sousa^^^®’^^^'’, 

C. Da Via*^, W. Dabrowski**^, A. Dafinca^^*^, T. Dai*^, O. Dale^^, E. Dallaire^^, C. Dallapiccola*®, 
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M. Dam^^, J.R. Dandoy^^, N.P. Dang^^, A.C. Daniells^^, M. Danninger^^^, M. Dano Hoffmann^^^, 

V. Dao^^, G. Darbo^*’®, S. Darmora^, J. Dassoulas^, A. Dattagupta^\ W. Davey^^, C. David'^^, 

T. Davidek^^^, E. Davies^^°’^ M. Davies^^^, P. Davison^*, Y. Davygora^^®, E. Dawe^*, L Dawson^^^, 
R.K. Daya-Ishmukhametova*^, K. De*, R. de Asmundis^°^®, S. De Castro^'^®’^*^’’, S. De Cecco***, 

N. De Groot^°^, P de Jong^*^’, H. De la Torre^^, E. De Eorenzi^^, E. De Nooij^°^, D. De Pedis^^^®, 

A. De Salvo^^^^, U. De Sanctis^^^, A. De Santo^^^, J.B. De Vivie De Regie^^’, WJ. Dearnaley^^, 

R. Debbe^^, C. Debenedetti^^^, D.V. Dedovieh®^, I. Deigaard^*^’, J. Del Peso^^, T. Del Prete'^^^’^^^*’, 

D. Delgove^^^, E. Deliot^^®, C.M. Delitzsch^^, M. Deliyergiyev^^, A. DeU’Aequa^*’, E. Dell’Asta^^, 

M. DeirOrso*^^^’^^^*’, M. Della Pietra^^^®’-^, D. della Volpe^^, M. Delmastro^, PA. Delsart^^, 

C. Deluea^**’, D.A. DeMareo*^*, S. Demers^’^, M. Demiehev^^, A. Demilly^*’, S.P Denisov^^*^, 

D. Derendarz^^, J.E. Derkaoui^^^‘^, E. Derue^°, P Dervan’^, K. Desch^\ C. Deterre^^, PO. Deviveiros^*^, 

A. Dewhurst^^^, S. Dhaliwal^^, A. Di Ciaccio^^^®’^^^'’, E. Di Ciaecio^, A. Di Domenico^^^®’'^^'’, 

C. Di Donato^°'‘^’i°^'’, A. Di Girolamo^o, B. Di Girolamo^o, A. Di Mattia^^^ B. Di Mieeoi 34 ad 34 b^ 

R. Di Nardo"*^, A. Di Simone^^, R. Di Sipio^^^, D. Di Valentino^^, C. Diaeonu*^, M. Diamond^^^, 

E. A. Dias^^, M.A. Diaz^^^, E.B. Diehl^^, J. Dietrich^^, S. Diglio*^, A. Dimitrievska^^, J. Dingfelder^\ 

P Dita^^®, S. Dita^^®, E. Dittus^*^, E. Djama*^, T. Djobava^^’’, J.I. Djuvsland^^®, M.A.B. do Vale^^‘^, 

D. Dobos^°, M. Dobre^®^, C. Doglioni^\ T. Dohmae^^^, J. Dolejsi^^^, Z. Dolezal^^^, 

B. A. Dolgoshein^^’*, M. Donadelli^^‘^, S. Donati^^^^’^^^*’, P Dondero^^^®’^^^'’, J. Donini^^, J. Dopke^^\ 

A. Doria^°^®, M.T. Dova’^, A.T. Doyle^^, E. Drechsler^^, M. Dris^*^, E. Dubreuil^^, E. Duehovni^’^, 

G. Duekeeki°°, O.A. Ducu^^^’^^ D. Duda^°’, A. Dudarev^o, E. Duliot^^ E. Duguid^’, M. Duhrssen^'’, 
M. Dunford^^®, H. Duran Yildiz^^, M. Duren^^, A. Durglishvili^^’’, D. Duschinger^^, M. Dyndal^^^, 

C. Eekardt^^, K.M. Ecker^°^, R.C. Edgar^^, W. Edson^, N.C. Edwards^^, W. Ehrenfeld^\ T. Eifert^®, 

G. Eigen*^, K. Einsweiler^^, T. Ekelof^^^, M. El Kaeimi^^^^^, M. Ellert^^^, S. Elles^, E. Ellinghaus^’^, 
A.A. Elliot^^^, N. Ellis^°, J. Elmsheuser^®*^, M. Elsing^'^, D. Emeliyanov^^^, Y Enari^^^, O.C. Endner^^, 
M. Endo^^^, J. Erdmann^^, A. Ereditato^^, G. Ernis'^^, J. Ernst^, M. Ernst^^, S. Errede^^^, E. Ertel*^, 

M. Esealier^^^, H. Esch^^, C. Eseobar^^^, B. Esposito^’, A.I. Etienvre^^^, E. Etzion^^^, H. Evans^\ 

A. Ezhilov^^^, E. Eabbri^*^®’^*^*’, G. Eaeini^\ R.M. Eakhrutdinov^^*^, S. Ealciano^^^^, R.J. Ealla’^, 

J. Ealtova^^^, Y Eang^^^, M. Eanti^^®’^^*’, A. Earbin^, A. Earilla^^^®, T. Earooque^^, S. Earrell^^, 

S. M. Earrington^’**, P Earthouat^*’, E. Eassi^^^'^, P Eassnacht^*^, D. Eassouliotis^, M. Eaucei Giannelli^’, 
A. Eavareto^*^®’^'^’’, E. Eayard^^’, P Eederie^^^^, O.E. Eedin^^^’"*, W. Eedorko^^^, S. Eeigl^®, 

E. Eeligioni^^, C. Eeng^^‘^, E.J. Eeng^, H. Eeng*^, A.B. Eenyuk^^°, E. Eeremenga^, 

P Eernandez MartinezS. Eernandez Perez^°, J. Eerrando^^, A. Eerrari^^^, P Eerrari^®^, R. Eerrari^^^®, 

D. E. Eerreira de Eima^^, A. Eerrer^^’, D. Eerrere^^, C. Eerretti^^, A. Eerretto Parodi^°^’^°'’, 

M. Eiasearis^\ E. Eiedler*^, A. Eilipcic^^, M. Eilipuzzi^^, E. Eilthaut^**^, M. Eincke-Keeler^^^, 

K. D. Einelli^^*^, M.C.N. Eiolhais^^^®’^^®‘^, E. Eiorini^^^, A. Eiran^°, A. Eischer^, C. Eischer^^, 

J. Eischer^’^, W.C. Eisher^*’, E.A. Eitzgerald^^, N. Elasehel^^, I. Eleek^^\ P. Eleisehmann^^, 

S. Eleischmann^^^, G.T. Eleteher*^^, G. Eletcher^^, R.R.M. Eletcher^^^, T. Elick^’^, A. Ploderus^\ 

E. R. Elores Castillo^®^, M.J. Elowerdew^^^, A. Eormica^^^, A. Eorti^^, D. Eournier^^^, H. Eox’^, 

S. Eracchia^^, P. Prancavilla^°, M. Eranchini^'^®’^*^’’, D. Eraneis^*^, E. Eraneoni^^^, M. Eranklin^’, 

M. Erate^^^, M. Eraternali^^'^’^^^*’, D. Ereeborn^*, S.T. Ereneh^^, E. Eriedrieh^^, D. Eroidevaux^®, 

J.A. Erost^^*^, C. Eukunaga^^^, E. Eullana Torregrosa^^, B.G. Eulsom^^^, J. Euster^^^, C. Gabaldon^^, 

O. Gabizon^’^ A. Gabrielli20^’20^ A. Gabrielli'32aa32b^ 5 Gadatsch^o^ S. Gadomski^^ 

G. Gagliardi^o^’SOb^ p Oagnon^i, C. Galea^o^ B. Galhardoi26^,126c^ p j Qallas^^o^ B j Qallopi^i, 

P. Gallus^^^ G. Galster^^, K.K. Gan“i, J. Gao^^b.ss^ Y. Gao^^ YS. Gao^^^.c^ PM Q^j-ay Walls'^^ 

E. Garberson^ 26 ^ p Garefa^^’^ j p Garcia Navarro^^’^ M Garcia-Sciveres^^, R.W. Gardner^^, 

N. Garelli^^^, V. Garonne^'^, C. Gatti^’, A. Gaudiello^°^’^°'’, G. Gaudio^^^®, B. Gaur^^\ E. Gauthier^^, 

P. Gauzzii32a,i32b^ J p Gavrilenko'^^ C. Gay^^^ G. Gaycken^i, E.N. Gazis^o, P. Ge^^^, z. Gecse^^^ 
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C. N.P. Gee^^^, D.A.A. Geerts^°^, Ch. Geich-Gimbel^\ M.P. Geisler^^®, C. Gemme^°^, M.H. Genest^^, 

S. GentileM. George^^, S. George^^, D. Gerbaudo^^^, A. Gersbon^^^, S. Gbasemi^^^, 

H. Gbazlane^^^'’, B. Giaeobbe^°^, S. Giagu^^^^’^^^*’, V. Giangiobbe^^, P Giannetti^^^®’^^^'’, B. Gibbard^^, 
S.M. Gibson’^ M. Gilebriesei^, T.PS. Gillam^^ D. Gillberg^o, G. Gilles^^, D.M. Gingrieb^-^, 

N. Giokaris'^, M.P Giordanii^^^’^^^c^ pjy[ Giorgi20^ F.M. Giorgi^^, PF. Giraud^^e^ p Giromini^’, 

D. Giugni‘’l^ C. Giuliani^^ M. Giulini5*^ B.K. Gjelsten^^^ S. Gkaitatzis^^^, f Gkialas^^^, 

E. L. Gkougkousis^*^, L.K. Gladilin^^, C. Glasman^^, J. Glatzer^*^, PC.F. Glaysber^^, A. Glazov^^, 

M. Goblirseb-Kolb^*^^, J.R. Goddard^^, J. Godlewski^^, S. Goldfarb^^, T. Golling^^, D. Golubkov^^*’, 

A. GomesR. Gonfalo^^^^, J. Gonealves Pinto Firmino Da Costa^^®, L. Gonella^\ 

S. Gonzalez de la Hoz^^’, G. Gonzalez Parra^^, S. Gonzalez-Sevilla^^, L. Goossens^*’, 

PA. Gorbounov^^, H.A. Gordon^^, I. Gorelov^°^, B. Gorini^*^, E. Gorini^^®’’^’’, A. Gorisek’^, 

E. Gornieki^^, A.T. Gosbaw^^, C. Gdssling^^, M.I. Gostkin^^, D. Goujdami^^^^^, A.G. Goussiou^^*, 

N. Govender^^^’’, E. Gozani^^^, H.M.X. Grabas^^^, E. Graber^^, I. Grabowska-Bold^^^, 

P. Grafstrdm^^®’^'^’’, K-J. Grabn^^, J. Gramling^^, E. Gramstad^'^, S. Graneagnolo^^, V. Grassi^^*, 

V. Gratebev^^^, H.M. Gray^°, E. Graziani'^^^, Z.D. Greenwood^^’", K. Gregersen’^, I.M. Gregor^^, 

P. Grenier^^^, J. Griffiths^, A.A. Grillo^^’, K. Grimm^^, S. Grinstein^^’", Pb. Gris^^, J.-E. Grivaz^^’, 

J.P Grobs^^, A. Grobsjean^^, E. Gross^’^, J. Grosse-Knetter^^, G.C. Grossi’^, Z.J. Grout^^^, E. Guan^^'’, 

J. Guentber^^^, E. Gueseini^^, D. Guest^^^, O. Gueta^^^, E. Guido^^^’^'^'’, T. Guillemin^^^, S. Guindon^, 

U. Gul^^, C. Gumpert"*^, J. Guo^^®, Y. Guo^^'’, S. Gupta^^°, G. Gustavino*^^^’^^^*’, P. Gutierrez^^^, 

N.G. Gutierrez Ortiz^^, C. Gutsebow^^, C. Guyot^^^, C. Gwenlan^^°, C.B. Gwilliam^^, A. Haas^^°, 

C. Haber^^, H.K. Hadavand^, N. Haddad^^^®, P. Haefner^^, S. Hageboek^^, Z. Hajduk^^, 

H. Hakobyan^’^, M. Haleem^^, J. Haley^^^, D. Hall^^**, G. Halladjian^*’, G.D. Hallewell*^, 

K. Hamaeber^’^, P. Hamal^^^, K. Hamano^^^, M. Hamer^^, A. Hamilton*^^^, G.N. Hamity^^^‘^, 

PG. Hamnett^^, E. Han^^'’, K. Hanagaki^^’^, K. Hanawa^^^, M. Hanee^^, P. Hanke^^®, R. Hanna^^^, 

J.B. Hansen^^, J.D. Hansen^®, M.C. Hansen^\ PH. Hansen^^, K. Hara^^*^, A.S. Hard^^^, 

T. Harenberg^’^, E. Hariri^^’, S. Harkusba^^, R.D. Harrington^^, PE. Harrison^™, E. Hartjes^®^, 

M. Hasegawa^’, S. Hasegawa^®^, Y. Hasegawa^^*^, A. Hasib^^^, S. Hassani^^^, S. Haug^^, R. Hauser^®, 

E. Hauswald^^, M. Havranek^^^, C.M. Hawkes^^, R.J. Hawkings^®, A.D. Hawkins^^ T. Hayasbi^^°, 

D. Hayden^®, C.P Hays^^**, J.M. Hays^^, H.S. Hayward’^, S.J. Haywood^^\ S.J. Head'^, T. Heek^^, 

V. Hedberg*^^, E. Heelan*, S. Heim^^^, T. Heim^’^, B. Heinemann^^, E. Heinrieb^^*’, J. Hejbal^^’, 

E. Helary^^, S. Hellman^^®^’^^^*’, D. Hellmieb^^, C. Helsens^^, J. HendersonR.C.W. Henderson’^, 

Y. Heng^’^, C. Hengler^^, A. Henriebs^’^, A.M. Henriques Correia^®, S. Henrot-Versille^^^, 

G.H. Herbert^^, Y. Hernandez Jimenez^^^, R. Herrberg-Sebubert^^, G. Herten^*, R. Hertenberger^®*^, 

E. Hervas^°, G.G. Hesketb^*, N.P Hessey^°^, J.W. Hetberly^*^, R. Hiekling^^, E. Higon-Rodriguez^^’, 

E. Hill^^^, J.C. Hill^^, K.H. Hiller^^, S.J. Hillier^*, I. Hinebliffe^^, E. Hines^^^, R.R. Hinman^^, 

M. Hirose^^^, D. Hirsebbuebl^^^, J. Hobbs^^*, N. Hod^°^, M.C. Hodgkinson*^^, P. Hodgson^^^, 

A. Hoeeker^*’, M.R. Hoeferkamp^*^^, E. Hoenig**’®, M. Hoblfeld^^, D. Hobn^^, T.R. Holmes^^, 

M. Homann^^, T.M. Hong^^^, E. Hooft van Huysduynen^^°, W.H. Hopkins^^^, Y. Horii^**^, 

A.J. Horton^^^, J-Y. Hostaeby^^, S. Hou^^\ A. Hoummada^^^^, J. Howard^^'^, J. Howartb^^, 

M. Hrabovsky^^^, I. Hristova^^, J. Hrivnae^^^, T. Hryn’ova^, A. Hrynevieb^^, C. Hsu^^^‘^, PJ. Hsu^^^’^, 

S. -C. Hsu'^^, D. Hu^^, Q. Hu^^’’, X. Hu^^, Y. Huang^^, Z. Hubaeek^^^, E. Hubaut^^, E. Huegging^^, 

T. B. Huffman^^°, E.W. Hugbes^^, G. Hugbes^^, M. Hubtinen^°, T.A. Hiilsing^^, N. Huseynov^^’^, 

J. Huston^*^, J. Hutb^’, G. laeobueei^^, G. lakovidis^^, I. Ibragimov^^^, E. leonomidou-Eayard^^^, 

E. Ideal^’^, Z. Idrissi^^^‘^, P. Iengo^°, O. Igonkina^°^, T. lizawa^’^, Y. Ikegami^^, K. Ikematsu^^^, 

M. Ikeno^^, Y. Ilebenko^^’'', D. Iliadis^^^, N. Ilie'^^, T. Inee^°^, G. Introzzi^^^®’^^^'’, P. loannou^, 

M. lodiee^^^®, K. lordanidou^^, V. Ippolito^’, A. Irles Quiles^^^, C. Isaksson^^^, M. Isbino^*, 

M. Isbitsuka^^^, R. Isbmukbametov^^^, C. Issever^^°, S. Istin^^^, J.M. Iturbe Ponee^^, R. luppa^^^®’^^^'’, 
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J. Ivarsson^', W. Iwanski^^, H. Iwasaki^^, J.M. Izen^^ V. Izzo^'^^^, S. Jabbar^, B. Jackson^^^, 

M. Jackson^^, P. Jackson^ M.R. Jaekel^*^, V. Jain^, K. Jakobs^^, S. Jakobsen^*^, T. Jakoubek^^’, 

J. Jakubek^^^, D.O. Jamin^^^, D.K. Jana^^, E. Jansen’^, R. Jansky^^, J. Janssen^^, M. Janus^’°, 

G. Jarlskog*\ N. Javadov^^’*, T. Javurek^*, L. Jeanty^^, J. Jejelava^^®’*, G.-Y. Jeng^^°, D. Jennens^^, 

P. Jenni^^’^ J. Jentzsch'^^, C. Jeske^’*^, S. Jezequel^, H. J. Jia^^^, Y. Jiang^^’’, S. Jiggins’^, 

J. Jimenez Pena^^^, S. Jin^^®, A. Jinaru^^®, O. Jinnouehi^^^, M.D. Joergensen^^, P Johansson*^^, 

K. A. Johns’, K. Jon-And^^^®’^^^'’, G. Jones^’°, R.W.L. Jones’^, T.J. Jones’^, J. Jongmanns^^®, 

PM. JorgeK.D. Joshi^^, J. Jovieevie^^^®, X. Ju^’^, C.A. Jung^^, P Jussel^^, A. Juste Rozas^’’", 

M. Kaei^^’, A. Kaezmarska^^, M. Kado^^’, H. Kagan^^\ M. Kagan^^^, S.J. Kahn^^, E. Kajomovitz^^, 

C. W. Kalderon^’°, S. Kama^°, A. Kamenshehikov^^'^, N. Kanaya^^^, S. Kaneti’^, V.A. Kantserov^^, 

J. Kanzaki^^, B. Kaplan^'®, E.S. Kaplan^’^, A. Kapliy^\ D. Kar^^, K. Karakostas^°, A. Karamaoun^, 

N. Karastathis^°’^*”, M.J. Kareem^^, M. Karnevskiy^^, S.N. Karpov®^, Z.M. Karpova^^, K. Karthik^^°, 
V. Kartvelishvili”, A.N. Karyukhin^^®, E. Kashif^’^, R.D. Kass^^\ A. Kastanas^^, Y Kataoka^^^, 

A. Katre"*^, J. Katzy^’, K. Kawagoe’*^, T. Kawamoto^^^, G. Kawamura^^, S. Kazama^^^, 

V.E. Kazanin^*^^’^, M.Y Kazarinov^^, R. Keeler^^^, R. Kehoe^°, J.S. Keller^’, J.J. Kempster”, 

H. Keoshkerian*^, O. Kepka^”, B.P Kersevan’^, S. Kersten^’^, R.A. Keyes*’, E. Khalil-zada^^, 

H. Khandanyan^^^^’^^^*’, A. Khanov^^^, A.G. Kharlamov^®^’*^, T.J. Khoo^*, V. Khovanskiy^’, 

E. Khramov^^ J. Khubua^i'’’", H.Y Kim*, H. S.H. Kim^^°, Y. Kim *b N. Kimura^^'^, 

O. M. Kind^^, B.T. King’"*, M. King*^’, S.B. King*^*, J. Kirk***, A.E. Kiryunin****, T. Kishimoto^’, 

D. Kisielewska**^, E. Kiss"**, K. Kiuehi*^**, O. Kivernyk**^, E. Kladiva*"*"**’, M.H. Klein**, M. Klein’"*, 

U. Klein’"*, K. Kleinkneeht**, P. Klimek*"*®^’*"*^*’, A. Klimentov**, R. Klingenberg"**, J.A. Klinger**^, 

T. Klioutehnikova***, E.-E. Kluge**^ P. Kluit***’, S. Kluth*°*, J. Knapik***, E. Kneringer“, 

E. B.E.G. Knoops**, A. Knue**, A. Kobayashi***, D. Kobayashi**’, T. Kobayashi***, M. Kobel"*"*, 

M. Koeian*"**, P. Kodys**^, T. Koffas*^, E. Koffeman***’, E.A. Kogan****, S. Kohlmann*’*, Z. Kohout***, 
T. Kohriki^^, T. Koi*"**, H. Kolanoski*^, I. Koletsou*, A.A. Komar^^ *, Y. Komori***, T. Kondo®^, 

N. Kondrashova"**, K. Kbneke"**, A.G. Kbnig***^, T. Kono^^, R. Konoplieh****’", N. Konstantinidis’*, 

R. Kopeliansky***, S. Koperny**®, E. Kbpke**, A.K. Kopp"**, K. Koreyl*^, K. Kordas**"*, A. Korn’*, 
A.A. Korol***^’"^, I. Korolkov**, E.V. Korolkova**^, O. Kortner****, S. Kortner****, T. Kosek**^, 

V. V. Kostyukhin**, V.M. Kotov^*, A. Kotwal"**, A. Kourkoumeli-Charalampidi**"*, C. Kourkoumelis^, 

V. Kouskoura**, A. Koutsman**^^, R. Kowalewski*^^, T.Z. Kowalski**^, W. Kozaneeki**®, 

A.S. Kozhin****, V.A. Kramarenko^^, G. Kramberger’*, D. Krasnopevtsev^*, M.W. Krasny***, 

A. Krasznahorkay***, J.K. Kraus**, A. Kravehenko**, S. Kreiss****, M. Kretz**'^, J. Kretzschmar’"*, 

K. Kreutzfeldt**, P. Krieger***, K. Krizka**, K. Kroeninger"**, H. Kroha****, J. Kroll***, J. Kroseberg**, 
J. Krstie**, U. Kruehonak^*, H. Kriiger**, N. Krumnaek^"*, Z.V. Krumshteyn^*, A. Kruse*’*, 

M.C. Kruse"**, M. Kruskal**, T. Kubota**, H. Kueuk’*, S. Kuday"**’, S. Kuehn"**, A. Kugel**^ 

E. Kuger*’"*, A. Kuhl**’, T. Kuhl"**, V. Kukhtin^*, Y. Kulehitsky^*, S. Kuleshov***’, M. Kuna***®’****’, 

T. Kunigo®*, A. Kupeo**’, H. Kurashige®’, YA. Kuroehkin^*, V. Kus**’, E.S. Kuwertz*^^, M. Kuze**’, 
J. Kvita***, T. Kwan*®, D. Kyriazopoulos**^, A. Ea Rosa**’, J.E. Ea Rosa Navarro*"*‘*, 
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